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Abstract: In this paper, we have made the attempt to classify the integrand basis of all two- loop diagrams 



in pure four-dimension space-time. Our classification includes the topology of two-loop diagrams which 
determines the structure of denominators, and the set of numerators under different kinematic configura- 
. tions of external momenta by using Grobner basis method. In our study, the variety defined by setting all 
CN • propagators to on-shell has played an important role. We discuss the structure of variety and how it splits 
to various irreducible branches when external momenta at each corner of diagrams satisfy some special 
• i— i . kinematic conditions. This information is crucial to the numerical or analytical fitting of coefficients for 
^ ■ integrand basis in reduction process. 
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1. Introduction 

In the last few years we have seen tremendous progress in the computation of one-loop amplitude (see 
reports [1,2] for reference) by using Passrino-Veltman(PV) reduction method [3]. The newly developed 
reduction methods can be sorted into two categories: the reduction done in the integral level, such as the 
unitarity cut [4-8] and generalized unitarity cut method [9, 10] and the reduction done in the integrand 
level, which was initiated by Ossola-Papadopoulos-Pittau(OPP) in [11] and further generalized in [12-17]. 
Comparing these two categories, the first one only focus on physical meaningful coefficients while the 
second one must also include spurious coefficients. Although only dealing with the integrand of loop 
integral, method of the second category is still very useful because we can easily program it and compute 
amplitude pure algebraically. 

Encouraged by the successful computation at one-loop level, it is natural to generalize these methods to 
higher loops, partially by our theoretical curiosity and partially by the requirement of prediction precision 
of cross section for experiments carried out on modern colliders such as Large Hadron Collider (LHC). 
However, things are much more involved comparing to one-loop case. The first difficult we are facing 
is that in general we do not know much about the basis for multi-loop amplitude. In fact, now it is 
clear that we should distinguish integral basis and integrand basis, and unlike the one-loop amplitude, the 
number of integrand basis is much larger than the number of integral basis in multi-loop amplitude. It is 
highly desirable to reduce integrand basis to integral basis further. One standard method of doing so is 
the Integrate-by-Part (IBP) method [18]. This method can be carried out in reasonable short time if the 
considered amplitude has only a few external particles, but it becomes unpractical with time consuming 
when the number of external particles increases. The second difficult we are confronting is how to extract 
coefficients of basis after knowing them. In other words, how can we perform the reduction for multi- 
loop amplitude practically? The IBP-reduction method will combine the task of finding basis and finding 
coefficients at the same time, while experience from one-loop amplitude calculation tells us that we can 
separate these two tasks. 

These difficulties in multi-loop amplitude computation have been addressed in the last few years by 
several works [19-29], where the main focus is on the reduction method at integrand level, which includes 
finding integrand basis and matching their coefficients. Initial step along this direction was taken in [19] 
where four-dimensional constructive algorithm for integrand has been applied to planar and non-planar 
contributions of four and five-point Maximal-Helicity-Violating(MHV) amplitudes in N = 4 Super- Yang- 
Mills theory. Using constraints from Gram matrix, similar determinant of monomials of numerators was 
achieved in [20]. Besides reduction in integrand level, reduction at the integral level is discussed in [21-25], 
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where in order to determine the physical contour for integral basis, the variety defined by setting all 
propagators on-shell has been carefully analyzed. 

Among these new developments, the application of computational algebraic geometry method to multi- 
loop amplitude calculation is very intriguing [27,28], where the Grobner basis plays a central role. It is 
quite easy to determine integrand basis by using Grobner basis method, although different linearly related 
set of integrand basis can be obtained with different choices of variable ordering in polynomial division. The 
coefficients of integrand basis can be determined by algebraic geometry method as well, and the knowledge 
of variety, including its branch structure and intersection pattern of branches, is very important in the 
application of this method. This algebraic geometry method has been tested by several examples in [27,29] 
at two and three-loop 1 . Encouraged by the success, in this paper we will use algebraic geometry method to 
systematically study all possible topologies of two-loop diagrams in pure four-dimension for any external 
momentum configurations, not only restrict to double-box or penta-triangle studied in various references. 
In our study, we will present the integrand basis and the structure of variety for all two-loop topologies. 

The paper is organized as follows. In section 2 we classify all possible topologies for two-loop diagrams. 
In section 3 the one-loop topologies are re-examined using algebraic geometry method, and some methods 
that will be used to analyze two-loop topologies are illustrated in these simple one-loop examples. In 
section 4, as a warm-up, we present the results of some trivial two- loop topologies. In section 5, a careful 
analysis of planar penta-triangle topology has been given, while in section 6, we give a detailed study 
of non-planar crossed double-triangle topology. In section 7, we summarize the results of all remaining 
topologies. In the last section, conclusions and discussions are provided. In appendix, we introduce some 
mathematical facts that can be used to study the branch structure of variety. 

2. An overview of general two-loop topologies 

In this section, we give an overview of general two-loop basis. Much of the results are scattered in literatures, 
and we assemble them here to make the paper self-contained. 

2.1 The general topology of two-loop diagrams 

Since two-loop diagrams can always be reduced to one-loop diagrams by cutting an inner propagator, we 
can inversely reconstruct two-loop ones by sewing two external legs of one-loop diagrams. The topology of 
one-loop diagram is very simple: we just attach various tree structures along the loop at some vertices Vi 
(see Figure [j] ) . 

From this diagram we can reconstruct two-loop diagrams by connecting two external legs, and there 
are several ways of doing so, which gives different two-loop topologies: 

• (A) If the two to-be-connected external legs are attached to the same tree structure, we will get 
two- loop diagrams as drawn in Figure ^. Explicit illustration shows that there are two connections. 

lr The numeric algebraic geometry method [30] can also be used if we only want to know the number of irreducible compo- 
nents. 
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Figure 1: The general topology of one-loop diagrams where various tree structures are attached along the loop at 
some vertices. 




(Al) (A2) 

Figure 2: The general topology of two- loop diagrams generated from one- loop diagrams by connecting two external 
legs attached to the same tree structure. The connection has been denoted by red color thick line. In connection 
(Al), two one-loop sub-structures do not share the same vertex while in connection (A2), they do share a common 
vertex. 

In the first connection (Al), two one-loop sub-structures do not share the same vertex while in the 
second connection (A2), they do share a common vertex. 

• (B) If the two to-be-connected external legs are attached to two nearby vertices along the loop, we will 
get two-loop diagrams as drawn in (B) of Figure |3[ All two-loop planar topologies can be generated 
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Figure 3: The general topology of two-loop diagrams of case (B) and case (C) by connecting two external legs 
attached to two different tree structures. For case (B), two tree structures are adjacent while for case (C), two tree 
structures are not nearby. 



from these diagrams. 

• (C) If the two to-be-connected external legs are attached to two non-nearby tree structures along the 
loop, we will get two-loop diagrams as drawn in (C) of Figure ||. All two-loop non-planar topologies 
can be generated from these diagrams. 

2.2 Classification of denominators of two-loop basis 

Having understood the general topology of two- loop diagrams, the next step is to classify the basis that can 
be used to expand any two-loop amplitudes. This is similar to the classification of scalar basis for one-loop 
diagrams, which includes box, triangle, bubble and tadpole. However, it is necessary to distinguish the 
integrand basis and the integral basis when talking about multi-loop amplitude basis. The integrand 
basis is that used by OPP method to expand expressions at integrand level, and some integrand bases will 
vanish after carrying out the integration. The integral basis is also called the master integrals(MIs). 
The number of integral basis is much smaller than the number of integrand basis, since after integration, 
some independent integrand bases can be linearly related to other integrand bases. The difference between 
these two kinds of basis can be easily seen in a simple example, the one-loop box topology where there is 
only one master integral fllfl2 j fl< , but two integrand bases: d 1 dIdzD 4 , an< ^ ^D^D^iof • Numerator of 
e • t with odd power will vanish after integration because of parity. 

To find the integrand and integral basis, we can use the procedure called PV-reduction. It involves 
many different manipulations on propagators and integral. Some manipulations are done at the level of 
integrand expressions, such as rewriting 2K\ ■£ = — (£— K\) 2 + £ 2 + K 2 . Some manipulations are carried out 
using the properties of integral, such as IBP method. Pure algebraic manipulations at the level of integrand 
expressions will produce the integrand basis, while combined with methods such as IBP, integrand basis 
can be reduced further to integral basis. 
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Above reduction has been discussed in many references and here we summarize some results (see, for 
example, [25] for details and reference) . For two-loop diagrams, denominators of integrand can always be 
written as products of three kinds of propagators 

V = DDD , (2.1) 

where 

D = lj(h - K ail ) 2 (h - K a rf...{t x - K a .n x -xf , 
D = £ 2 (£ 2 - K bjl ) 2 (£ 2 - K by2 ) 2 ...{£ 2 - if M2 „i) 2 , 

D = {£ 1 + £ 2 + K^) 2 ...^ + £ 2 + K C:ns ) 2 , (2.2) 

and n±, n 2 are numbers of propagators that containing only l\ or £ 2 , while n 3 is the number of propagators 
containing £\,£ 2 at the same time. By the symmetry among three kinds of propagators (i.e., the relabeling 
of momenta £\,£ 2 ), we can always restrict the number of propagators with condition 

n\ > n 2 > n 3 • (2-3) 

The up-bound of ni,n 2 ,n^ and their summation depends on the space-time dimension of problem. For 
example, if we consider physics in (4 — 2e)-dimension, we could have 

ni, n 2 , n 3 < 5 , n\ + n 2 + n 3 < 11 . (2.4) 
But if we constraint to pure four-dimension, the condition becomes 

ni, n 2 , n 3 < 4 , rii + n 2 + n 3 < 8 . (2.5) 
It is possible to classify denominators of integrand and integral basis by combining conditions ( |2.3[ ) and 



(2.5) for 4-dimension case (or ( |2.3D and fl2.4| ) for (4— 2e)-dimension case). For (4 — 2e)-dimension, conditions 



(|2.3|) and ( fO| ) requires n 3 < 3. Thus if we arrange all possible solutions of (ni,ri2,n 3 ) by value of n 3 , we 



have following 4 groups of solutions 

n 3 = 3: (5, 3, 3), (4, 4, 3), (4, 3, 3), (3, 3, 3) ; 

n 3 = 2 : (5, 4, 2), (5, 3, 2), (4, 4, 2), (5, 2, 2), (4, 3, 2), (4, 2, 2), (3, 3, 2), (3, 2, 2), (2, 2, 2) ; 
n 3 = 1 : (5, 5, 1), (5, 4, 1), (5, 3, 1), (4, 4, 1), (5, 2, 1), (4, 3, 1), (5, 1, 1), (4, 2, 1), 

(3, 3, 1), (4, 1, 1), (3, 2, 1), (3, 1, 1), (2, 2, 1), (2, 1, 1), (1, 1, 1) ; 
n 3 = : (5, 5, 0), (5, 4, 0), (5, 3, 0), (4, 4, 0), (5, 2, 0), (4, 3, 0), (5, 1,0), (4, 2, 0), (3, 3, 0), 

(4,1,0), (3,2,0), (3,1,0), (2,2,0), (2,1,0), (1,1,0) . (2.6) 
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For pure four-dimension, the number of solutions decreases a lot, since instead we have 713 < 2, and the 
possible solutions of (m, n 2 ,ri3) for ([T^) and (|2.5D are listed into 3 groups: 



n 3 = 2 : (4, 2, 2), (3, 3, 2), (3, 2, 2), (2, 2, 2) ; (2.7) 

n 3 = 1 : (4, 3, 1), (4, 2, 1), (3, 3, 1), (4, 1, 1), (3, 2, 1), (3, 1, 1), (2, 2, 1), (2, 1, 1), (1, 1, 1) ; 

n 3 = : (4, 4, 0), (4, 3, 0), (4, 2, 0), (3, 3, 0), (4, 1, 0), (3, 2, 0), (3, 1, 0), (2, 2, 0), (2, 1,0), (1, 1, 0) . 

Topologies with 71,3 = contain two-loop diagrams coming from sewing two one-loop structures at one single 
vertex as shown in Figure || , while topologies with 713 = 1 contain planar two-loop diagrams with one 
common propagator shared by the two one-loop structures as shown in Figure ||. All two-loop non-planar 
diagrams are included in 77,3 = 2 group as shown in Figure ||. 

The topologies of two-loop basis can be distinguished by (774, 77,2, 723), while in order to get the necessary 
integrand or integral basis, we should also dress propagators with numerators. In two-loop case, the so 
called scalar integrand(with numerator 1) is not enough to expand all diagrams, we need also terms with 
numerators that are constructed from monomial of scalar products of external and internal momenta. The 
distinction between integrand and integral basis becomes important when discussing the classification of 
numerators, and in this paper, we will focus only on the integrand basis in pure four-dimension. 

3. The integrand basis of one-loop diagrams in pure four-dimension 

As a warm-up, we take the integrand basis of one- loop diagrams as a simple example to demonstrate 
various ideas that we will meet in later part of this paper. All results in this section are known in other 
references such as [11,27,28], however, we recall them here for self-contained since these results are also 
related to the integrand basis of two- loop diagrams with 77,3 = 0. 

In pure four-dimensional case, since each external or internal momentum has four degrees of freedom, 
we need four independent momenta to expand all kinematics. One way of constructing the momentum 
basis is from two arbitrary independent momenta K\,K 2 . Assuming (K\ + K 2 ) 2 7^ 0, we can construct 
following four null momenta Ci, i = 1, 2, 3, 4 



■ , „ K\ + K 1 -K 2 - sgn(K 1 • K 2 )\ y^g • K 2 f - K'{K'j\ 

61 = _ 1 Kl Wx + K 2 y Kl2 




» , , Kl + K l -K 2 - sgn(K 1 • K 2 )\^{K X ■ K 2 f - ggg \ 

62 = _ 1 K2 (K 1 + K 2 y Kl2 j ' 

63 = — 2i~ ' 64 = — a" ' (3 ' 1} 

where ^\ 2 = ^^(j^^jf^ 1 ^" 2 ^ ■ This momentum basis has following property: among all inner products of 
e, • Bj, the only non-zero ones are e\ ■ e 2 = 1 and ■ = l 2 . Definition ( |3.1[ ) also makes massless limit 

2 In fact, this property has not determined d uniquely, since there is a freedom to redefine e\ —¥ wet and e2 — > w~ 1 e2 and 
similarly for e3,e4 pair. 
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smoothly, i.e., when K\ — > 0, e\ ~ i^i and when — )• 0, e 2 ~ if 2 . Using above momentum basis, we can 
expand any external and internal momenta, such as 

4 

Ki = (Ki • e 2 )ei + (Ki • ei)e 2 + (ifj • e 4 )e 3 + (iQ • e 3 )e 4 = ^ a^e, , 

£={£■ e 2 )ei + • ei)e 2 + • e 4 )e 3 + • e 3 )e 4 = x 2 ei + xie 2 + x 4 e 3 + x 3 e 4 , (3.2) 
and the Lorentz invariant scalar products are given by 

4 

I 2 = X\X 2 + X 3 X 4 , i ■ Ki = a ij x j ■ (3-3) 

The importance of above expansions ([T^) and ( |3.3[ ) is that any integrand expression can be written 
as a rational function rt ^n)f 2 f 3 f 4 \ 1 1 an d the PV-reduction procedure is equivalent to finding following 
expansion of numerator 

/( x .) = ^ A(^) + r(xi) , (3.4) 
t 

where the remaining polynomial r(xi) is nothing but the integrand basis we are looking for. In a more 
mathematical language, propagators D t generate an ideal I in polynomial ring k[x%, x 2 , x 3 , x 4 ], and the 
integrand basis is constructed by representative elements in the quotient ring k/I under some physical 
constraints. One physical constraint we will use is the total degree ri£ of loop momentum I in numerator. 
For renormalizable theory, we have n,£ < rip where no is the number of propagators in denominator. 

Having these general preparations, we move to explicit discussion of various one-loop topologies, such 
as box, triangle, bubble and tadpole [11,27,28]. For simplicity, we will consider massless propagators only, 
but the massive ones can be discussed in similar way. 

3.1 One-loop box topology 

For box topology, four propagators are given by 

D = £ 2 , D 1 = (£- K x f ,D 2 = {£-K 1 - K 2 f , D 3 = (£ - K t - K 2 - K 3 ) 2 . (3.5) 

Without loss of generality, we can use K\ , K 2 to construct a set of momentum basis and expand all 
momenta. There are 4 variables (x\, x 2 , x 3 , x 4 ) coming from loop momentum £, and all above propagators 
can be translated into following polynomials of Xi variables 

D = x\x 2 + x 3 x 4 , Di = D - 2(anxi + a i2 x 2 ) + 2a\\a\ 2 , 
D 2 = D — 2(a 2i xi + a 22 x 2 ) + 2a 21 a 22 , 

D 3 = D - 2(a 3 ixi + a 32 x 2 + a 33 x 3 + a 34 x 4 ) + 2a 3i a 32 + 2a 33 a 34 , (3.6) 
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where we have used the parametrization K\ + ... + K \ = Y2t=l a n e t- It ^ s eas Y to see that (Dq — D\)/2 = 
oinxi + 012^2 — ca\\a.Yi belongs to the ideal generated by Do, D\, D2, D3. However, the linearity of this 
equation means that in quotient ring k[x±, X2, X3, x^j/I, we can always treat variable x\ as combination 
012 — ^77 x 2- In other words, we can use equation = a.\\x\ + 012X2 — 011012 to solve x\ and eliminate 
variable x\ from the quotient ring k[xi, X2, x 3 , x^jl . Similarly we can use other two linear equation 
(D — D 2 ), (D — D 3 ) to solve, for example, variables X2,x% and get 

012022(011 - «2i) ana 2 i(ai2 - 022) 
x\ = , x 2 = , 

ailO:22 — "12021 "12021 — 011022 

-x 4 a 34 - «2i032 + 031032 . 012(011 - a2i)(a 2 ia 32 - 022031) 

X 3 = 1 r h Q34 • (3.7) 

033 033(Oll«22 - Oi 2 02l) 

Since x\, x 2 , x 3 can be solved as linear polynomial of X4, we will call them reducible scalar products (RSP), 
while the remaining variable X4 will be called irreducible scalar products (ISP). 

After substituting solution ( p.7| ) into Dq we get a quadratic polynomial of variable X4 

D (x4) = — — —x\ + C1X4 + c , (3.8) 
033 

where 

012021(031(032 - Q22) + 0330:34) + 011(012(022031 - 021032) + o 2 2((o2i - 031)032 - 0330-34)) 



ci 



(012021 - 011022)033 



011012021022(011 - a 2 l)(oi2 - 022) ,„ n , 

c = • (3.9) 

(012021 - 011022) 



The problem is then reduced to quotient ring k[x^\/ (D^x^)). Since ( |3.8[ ) is a quadratic polynomial, the 
representative elements in quotient ring can take following two terms: 1 and X4. It is worth to notice that 
although in this example the quotient ring is finite, it is not true in general. In fact, if we consider the 
quotient ring as linear space, in general the dimension of it will be infinity, i.e., there are infinite number 
of representative elements. Only when some constraints are taken we get finite number of representative 
elements. 

There is another issue regarding to the ideal defined by ( |3.8|) . The quadratic polynomial is reducible, 
i.e., it can be factorized as product of two factors a(x4 — zi)(x4 — 22) where z%,Z2 are two roots. This 
will split the solution space into two branches, which are obtained by setting either factor to zero. The 
variety 3 defined by this polynomial is the union of two branches (here is just two points). Both branches 
are needed to analytically (or numerically) determine coefficients of two integrand basis £) Q DiD2D 3 an( ^ 
r) £^ A £, 2 Dz a ^ ^ e integrand level. One of the main focus of this paper is the study of variety determined 
by setting all propagators of a given topology to zero, and their branch structures as well as degeneracy for 
special configurations of kinematics, such as massless limit of external momenta or some attached momenta 
becoming zero. 



3 we call the solution space as variety, which is the terminology used in algebraic geometry. 
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3.2 One-loop triangle topology 



The three propagators are given by Dq, D\, D2 as in (3.5), thus we can solve 

ai2«22(an - «2i) "11021(022 - "12) ,„.„> 
x\ = , x 2 = . (3.10) 

"11^22 — 012021 CH1022 — «12«21 

For triangle, x%, X2 become RSPs, and only X3,X4 are left as independent variables. Putting them back to 
Do we obtain 

n , \ anai2a2ia22(an - a2i)(«22 - "12) . /oun 

^0(^3,^4) = 7 y> r-x 3 x 4 . (3.11) 

(012021 - an«22 ) 

The quotient ring is given by k[x%, X4]/ {Dq(x^, X4)), and its representative elements can be taken as 
l,Xg 3 ,X4 4 with 713,714 > l 4 . Unlike the box topology, the dimension of this quotient ring will be infinity. 
To select finite number of representative elements from quotient ring, we constraint the power 723, 77,4 of 
X3 3 , X4 4 no more than three. This corresponds to the condition that the power of £ in numerator is 
no more than three for triangle topology. In this case we get following seven representative elements 
{1, X3, X4, x|, £4, £3, X4} which were given in [11]. 

After getting the integrand basis, we need to consider the problem of fitting their coefficients. For 



this purpose, understanding the variety defined by ( 3.11 ) becomes important. Assuming that the equation 
is given by X3X4 + d = with d ^ 0, we can solve X3 = — d/x 4 . Putting it back to integrand basis we 
can rewrite them as seven monomials of X4 only: x\ with t = —3,-2,-1,0,1,2,3. Thus to find their 
coefficients, we just need to substitute xi,X2,X3 as functions of X4 into integrand expressions obtained by 
Feynman diagrams or sewing three on-shell amplitudes using unitarity cut method. Having the expressions 
of X4, we can identify corresponding coefficients for a given power of X4. For numerical analysis, we can take 
seven arbitrary points of X4 to construct seven linear equations and find the seven unknown coefficients of 
integrand basis. 

There is a technical issue regarding the method we have just described. To guarantee that we will 
get exactly the form Ylt=-3 c t x % we mus t first subtract all contributions from box topologies. Similar 
manipulation should be taken when fitting coefficients for bubble and tadpole at one-loop. In other words, 
we should subtract contributions from all other higher topologies which contain the same set of propagators 
in the problem. 

The procedure we have just described is called parametrization of variety. For the simple example with 
d ^ 0, there is only one irreducible branch parameterized by X4. However, with some special kinematic 
configurations, above irreducible branch can split to two branches. This happens when Kf = 0, so a\2 = 
or K% = 0, so («2i — an)(«22 — 012) = 0, or K% = (K\ + K2) 2 = 0, so 021^22 = 0. In other words, 
when at least one leg is massless, the definition equation of variety is reduced to X3X4 = 0, and we get two 
irreducible branches. The first branch is parameterized by setting X3 = with X4 as free parameter, and the 
second branch is parameterized by setting X4 = with X3 as free parameter. If we put the parametrization 



Using (3.11) we can eliminate any product of £3,2:4. 
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of the first branch back to integrand basis, x% with n = 1,2,3 will be zero and their coefficients can not 
be detected by method described in previous paragraph, which means that the first branch can only be 
used to find four coefficients of integrand basis. Similarly, the second branch can only be used to find 
coefficients of integrand basis 1, x 3 , x%, x|. These two branches intersect at one point X3 = X4 = 0, thus we 
have 4 + 4 — 1 = 7, i.e., both branches are necessary to fully determine coefficients of integrand basis and 
they provide the complete information of all coefficients. 

3.3 One-loop bubble topology 

There are two external momenta, and because of momentum conservation only one of them is independent. 
In this case, we take K\ and another auxiliary momentum P to construct the momentum basis. With two 
propagators Do > D\ we can solve 

xi = a 12 (l - —) , (3.12) 

and there are three ISPs (X2, x 3 , £4). After eliminating xi, Dq becomes polynomial of three ISPs 

Do{x 2 -, x 3 , 24) = x 2 au ( 1 - — ) + 23x4 , (3.13) 

V "11/ 

which defines the variety of ideal generated by Dq in polynomial ring k [x 2 , x 3 , X4] . Unlike box and triangle 
topologies, it is hard to find representative elements in the quotient ring 

k[x 2 ,x 3 ,x 4 ]/ {D (x 2 ,x 3 ,x 4 )) 

and we need a systematic way to do so. A good solution for this problem is to use the Grobner basis of 
ideal. Firstly we write down all possible monomials 4 4 with n 2 + n 3 + 724 < 2 to fulfill physical 

constraints, and then divide each monomial X 2 X ^ X ^ by Grobner basis and collect all monomials in the 
remainder. These monomials collected from the remainder times D ^ Dl give the integrand basis. 

A technical issue of above algorithm is the ordering of variables in the constructing of Grobner basis. 
Different ordering gives, in general, different sets for representative elements and they are equivalent to 
each other. Once a particular ordering is chosen, we should stick to it through the whole calculation to 
avoid inconsistency. For instance, if the ordering is chosen as X3 > X4 > x 2 we get 9 integrand basis as 

{1,X 2 ,X 2 2 ,X3,X 2 X 3 ,X3 2 ,X4,X 2 X4,X 4 2 } . (3.14) 

This integrand basis can be used to expand bubble topology. In order to get the coefficients of integrand 
basis analytically, we should first put t = x 2 e± + x\e 2 + x 4 e 3 + x^e^ back into integrand expression obtained 
from Feynamn diagrams or unitarity cut method after subtracting all box and triangle contributions. Then 
we can replace xi by ( |3.12 ), and get a polynomial A(x 2 ,X3,X4) with three ISPs (x 2 ,X3,x 4 ). The next step 



is to divide this polynomial by Gronber basis and obtain the remainder. This algorithm, different from 
previous parametrization method, ensures that the remainder is nothing but the linear combination of 
monomials in integrand basis with coefficients we want to find. 
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^5 „,„ _ ^2 n ^2 

as well as integrand basis, and get 



If using the parametrization method, we can replace 5 x% = — ai 2 f|(l— jjp-) in the expression A(x 2 , X3, X4) 



c(6)x 2 4 af 2 2c(6)x- 2 3 af 2 c(5)x 2 3 ai 2 c(6)x 2 2 af 2 c(4)rr 2 2 a 12 c(5)x 2 2 a i2 
x^a^ X4 2 an x^au X4 2 x^acu X4 

_ c(4)x 2 a 12 + C ( 3 ) X2 2 + C ( 8 ) X2X4 + c (2)a; 2 + c (9)x 4 2 + c(7)x 4 + c(l) . 
X4 

Since we have already used the equation to reduce one variable further, remaining variables x 2 ,X4 are free 
variables and there is no need to divide Grobner basis further. What we need to do is to compare each 
independent monomial x\x\ (a, b could be negative integers) at both sides. The parametrization can also 
be used for numerical fitting, and we only need to get enough linear equations to solve coefficients by taking 
sufficient numerical points (22 > £4) at both sides. 

Similar to triangle topology, the variety defined by ( |3.13 ) is irreducible for general kinematic config- 



uration. However, when K\ = 6 we have a,\i = by our construction, thus equation ( [3,13 ) is reduced 



to X3X4 = 0. In other words, the irreducible variety is degenerated to two branches: one with X3 = and 
x 2 ,X4 as free parameters, and another with X4 = and x 2 ,X3 as free parameters. Each branch can detect 
six coefficients out of nine integrand basis, while three basis {l,x 2 ,x|} can be detected by both branches. 
Thus we have 6 + 6 — 3 = 9, and both branches are necessary to find all coefficients of integrand basis 
analytically or numerically. 

3.4 One-loop tadpole topology 

In this case, we choose arbitrary two independent momenta to construct the momentum basis. Since there 
is only one propagator Dq, all four variables Xi, i = 1, 2, 3, 4 are ISPs and the variety is defined by equation 

D = x\ x 2 + £3X4 . (3.15) 

Requiring the total dimension of monomials no more than one, we get following basis 

{l,xi,x 2 ,x 3 ,x 4 } . (3.16) 

This variety is irreducible and we can parameterize it by solving x\ = —^ A . Thus after putting x\ back 
to expressions obtained from Feynman diagram after subtracted all contributions from boxes, triangles 
and bubbles, we can read out coefficients of one-loop tadpole integrand basis by comparing monomials of 
x 2 x 3 x 4 . 



4. A premiere: some trivial two-loop topologies 

Starting from this section, we will discuss the integrand basis and variety of various two-loop topologies 
classified in (|2.7p using the same method presented in previous section for one-loop diagrams. Before going 

5 This parametrization works for almost every value of £2 except X2 = and X2 = an where the variety is degenerate. 
6 For one-loop theory, bubble basis with K\ — vanishes after integration, but it is necessary at the integrand level. 
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(A31) (A22) (A21) (All) 

Figure 4: The type (A) contains all 10 diagrams with 713 — 0. Every diagram is denoted by (Anm) where n, m is 
the numbers of propagators of the left and right one-loop diagram respectively. These diagrams are drawn in most 
general form, and some external momenta, for instance if 7 in (A44) diagram, could be absent. All external momenta 
are out-going while convention of loop momenta is labeled as an arrow in each diagram. 

to non-trivial topologies, there are some topologies whose integrand basis and structure of variety are quite 
simple. These includes all topologies of type (A), of which two one-loop structures share only one single 
vertex, and the topologies that sharing common propagators between two one-loop structures but have 
maximal number of propagators, i.e., 8 propagators for pure 4-dimensional two-loop diagrams. 

4.1 Two- loop topologies of type (A): single-point attached diagrams 

All two-loop topologies of type (A) can be found in Figure (f|). Since there is no propagator involving 
at the same time, integrand basis and variety defined by propagators will be double copy of corresponding 
two one-loop sub-diagrams with minor modification. This modification comes from constraints of total 
degree of monomials in integrand basis. Taking topology (A33) as an example, we can use if 1, if 2 to 
construct momentum basis ei,e2, es,e^, and X3 = i\ ■ e^, X4 = £% ■ will become ISPs after solving 
linear equations. Similarly we can use if 3, if 4 to construct another momentum basis ei,e2,e3,e4, and 
2/3 = (-2 "63, 2/4 = £2 -&4 become ISPs. Thus the representative elements of (A33) can be given by monomial 
x z ix A i yT Z yT A with constraint 713 + 714 < 3 because along the loop there are three vertices. Similarly we 
have 7713 + 1714 < 3. However, since there are only five points along whole two-loop diagram we should have 
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n% + m + 7773 + ttt-4 < 5. With this condition should be excluded from integrand basis and we get 
7 x 7 - 4 = 45 basis. 

The variety is also the union of varieties of corresponding two one-loop sub-diagrams, so its structure 
can be easily inferred. To determine coefficients of integrand basis, similar procedures as presented in 
previous section can be applied, such as using Grobner basis method or parametrization method. 

4.2 Topologies with maximal eight propagators 

Besides topologies of type (A), there are three special topologies in type (B) and (C) which have maximal 
eight propagators. Since there are in total eight degrees of freedom for two- loop momenta £i,£2, on-shell 
conditions of eight propagators will completely freeze all eight freedoms, and the variety will be fixed to 
isolated points. These topologies are planar penta-box as shown in (B43) of Figure (||), and non-planar 
crossed penta-triangle, crossed double-box as shown in (C42), (C33) of Figure (H) 7 . Their results are pretty 
simple and we will present them in this subsection. 

4.2.1 The topology (B43): planar penta-box 

For (B43) topology, we take K\,K§ to construct momentum basis ej, i = 1,2,3,4 and expand both loop 
momenta £±,£2 with coefficients x\ = £\- e\ and y^ = £2 ■ e%- Since there are four propagators containing 
only £\, just like the one-loop box case, x\, X2, X3 can be solved as linear functions of X4. After substituting 
these solutions, Do = £\ becomes quadratic function of single variable £4 

Dq = c 2 xl + C1X4 + c , (4.1) 

where Cj are some functions of external momenta, which may be complicated according to the explicit 
kinematic configuration, but not important here. Similarly, there are three propagators containing only 
£2, so like the one-loop triangle case, 2/1,2/2 are solved as linear functions of 2/3,2/4. After substituting these 
solutions, Dq = £\ becomes 

A) = C202/I + C02Z/4 + CU2/32/4 + C10Z/3 + 0)12/4 + coo • (4.2) 

Propagator {£\ + £2 + Kq) 2 can also be expressed as function of these ISPs as 

Do = ^2 dijx\y 3 3 + ^2 dijx\y{ , ij = 0,1 , (4.3) 

ij ij 

where we have used the conditions £\ = £\ = 0. 

The integrand basis is constructed by dividing monomials x^y^y™ 4 with conditions 724 < 5, 711,3 + 
7714 < 4 and 714 + 7713 + 7114 < 7 over Grobner basis of the ideal generated by polynomials Do, Do, Do- The 
result is 

#B43 = {l,£4,2/3,2/4} • (4.4) 

7 Topology (A44) also has eight propagators. The variety is simply given by four isolated points and integrand basis has 
exactly four terms. These four points can be used to determine four coefficients of basis. 
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The variety defined by Dq,Dq,Dq has four branches and each branch has a single solution. Thus using 
four branches, we can fit coefficients of four integrand basis analytically or numerically by the method 
discussed in previous section. 

Above results will not change whether if Kq, Kj are absent or some of K{, i = 1, 2, 3, 4, 5 are massless. 

4.2.2 The topology (C42): non-planar crossed penta-triangle 

For (C42) topology, we take K\,K± to construct momentum basis ei,i = 1,2,3,4 and expand both loop 
momenta i\ , £2 into this basis with coefficients X{ = i\- e« and yi = £2 ■ e« . Since there are four propagators 
containing only £\, x\,X2,x% can be solved as linear functions of X4, and Dq = £\ can be rewritten as a 
polynomial of X4 

Dq = C2X\ + C1X4 + c . (4.5) 

Coefficients Cj are again some functions of external momenta whose explicit expressions are not important 
here. Similarly, there are two propagators containing only £2, and 2/2 can be solved as linear function of 
2/1,2/3,2/4- However, unlike the topologies of type (B), here we have two propagators containing both £±,£2, 
i.e., D = {£\+£2 + Kq) 2 and D\ = (£i + £2 + K 6 + K 7 ) 2 . We can get one more linear equation D\ — Dq and 
solve yi as linear function of £4, 2/3, 2/4. Thus we have three ISPs (0:4, 2/3,2/4) and three quadratic polynomials 
with respect to these three ISPs. Using ideal generated by three polynomials we find the integrand basis 
is given by 

&C42 = {1,^4,2/3,2/4} . (4.6) 

The variety defined by these three quadratic equations has four branches, and each branch is given by a 
point. Thus using four branches we can find coefficients of four integrand basis. Again above discussion 
does not change whether if Kq,K 7 are absent, or any of other external momenta go to massless limit. 

4.2.3 The topology (C33): non-planar crossed double-box 

For (C33) topology, we take K\,K± to construct momentum basis e«,i = 1,2,3,4, and expand both loop 
momenta £\,£2 into this basis. It is possible to get five linear equations from the eight on-shell cut equations, 
and solve, for instance, xi, X2, X3, 2/1, 2/2 as functions of three ISPs (24, 2/3, 2/4). After substituting all RSPs 
in the remaining three propagators we get three quadratic polynomials with respect to three ISPs. The 
variety defined by these three quadratic polynomials is given by eight points, and by using Grobner basis 
method, we can find integrand basis with 8 elements 

#C33 = {1, X 4 , 2/3, xm, 2/3, 2/3, 2/4, 2/32/4} • (4.7) 

As usual, each point of variety can detect one coefficient of integrand basis, and by using all 8 points, 
we can get all coefficients. Again above discussion does not rely on the explicit configuration of external 
momenta. 
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(B32) (B31) (B22) (B21) (Bll) 

Figure 5: The type (B) contains 9 diagrams with 77,3 = 1. Every diagram is denoted by (Bnm) where n, m are the 
numbers of propagators containing pure t\ or £2 loop momentum. The diagrams are drawn in most general form, 
and some external momenta, for instance Kq^Ki in (B43), could be absent. All external momenta are out-going 
while convention of loop momenta is labeled as an arrow in each diagram. 

5. Example one: planar penta-triangle 

Having understood simple topologies of planar penta-box, non-planar crossed penta-triangle and crossed 
double-box, let us move to non-trivial topologies where varieties are given by manifold with dimension at 
least one, not just isolated points. For these topologies, analysis becomes more complicated, so we will 
take two topologies as examples to illustrate various properties. The first example we will study is planar 
two-loop penta-triangle topology, as shown in (B42) of Figure [|. 

The penta-triangle topology has 7 propagators. If we choose K\,K& to generate momentum basis 
ei, i = 1, 2, 3, 4, all kinematics can be expanded as 

4 = x 2 ei + xie 2 + x 4 e 3 + x 3 e 4 , l 2 = y 2 e.\ + y±e 2 + y 4 e 3 + 7/ 3 e 4 , 
K\ = a\\e\ + a\ 2 e 2 , K\ + K 2 = a 2 \e\ + 02262 + «23<3 3 + «24e 4 , 
K\ + K 2 + K 3 = a 3i ei + a 32 e 2 + a 33 e 3 + a 34 e 4 , 

K A = /3nei + (3i 2 e 2 , K 5 = jud + 7i 2 e 2 + 7i 3 e 3 + 7i 4 e 4 , (5.1) 

and Kq is constructed from momenta conservation. Above parameters are general if Kq is arbitrary, but 
when Kq = or Kq is massless, there will be relations among parameters. For example, when Kq = 0, we 
should have 

711 = -fill ~ "31 , 712 = ~/3l2 - "32 , 713 = -"33 , 714 = ~«34 , (5.2) 

and if K$ = Kq = at the same time we should have jn = and 

"31 = -/?n , "32 = -A2 , «33 = , a 34 = . (5.3) 
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These relations will be important when discussing branch structure of variety under special kinematic 
configurations. 

Using above expansion, we can expand all the following seven propagators 

D Q = l\, D^ih-K^ 2 , D 2 = {£ 1 -K 1 -K 2 ) 2 , D 3 = {£ 1 -K 1 -K 2 -K 3 ) 2 

D = £l D 1 = (£ 2 -K 4 ) 2 , D = (£ 1 + e 2 + K 5 ) 2 , (5.4) 

and use four linear equations D\ — Dq = 0, D 2 — Do = 0, D 3 — Do = and D\ — Dq = to solve x±, x 2 , x 3 , y 2 
as linear functions of four ISPs (X4, y±, y 3 , 7/4), which are given by 

_ Ql2("24«33 - Q23Q34) 

"12("23«31 - 021^33) + aii(a 2 2"33 ~ 23 a 32 ) 



(5.5) 



| Ql2(-Q21«22Q33 + «ll(«22«33 - 023032) + «23(«31«32 + Q33Q34 ~ Q33Q24)) 

ai 2 (a23"3i - "21033) + an(a22«33 - "23032) 

_ «ll(«23Q34 - «24«33) 

ai 2 (a 2 3 a 3i — "21033) + 011(022033 - o 23 q 32 ) 

l 011(021022033 + 012(023031 - 021033) - 023(031032 + 033034 - 033024)) ^ 5 ^ 

oi2(o 23 a 3 i - a 2 ia 33 ) + an(a 22 a 3 3 - a 23 a 32 ) 

_ 012(021034 - 024031) + 011(024032 - 022034) 
oi2(o 23 a 3 i - a 2 ia 33 ) + an(a 22 a 33 - 023032) 

| Oll((-023024 - 022021 + 02203l)032 + Ol2(o2lQ32 ~ O22O31) + O22O33O34) 

oi2(a 23 a 3 i - a 2 ia 33 ) + aii(a 22 a 3 3 - a 23 a 32 ) 

| 012(023024031 + 021(031022 ~ 031032 ~ 033034)) ^ ^ 

oi 2 (a 2 3«3i - a 2 i«33) + on(Q 22 a 3 3 - a 23 a 32 ) 



and 



y 2 = /3n(l-| L ) • (5.8) 
P12 

Now we consider the remaining three equations. Firstly the equation Dq = becomes a quadratic 

r r 

equation of X4 and we always have 2 solutions x 4 1 ,x 4 2 in C-plane. There are no intersections between 
these two solutions, so the variety defined by this equation has been split into two separate branches 
parameterized by x\. Remaining two equations are 

Do = fti (1 - | L )yi + ym = , (5.9) 

Pl2 

,r i / r , ->Pii\„. , / r , \„. , / r 



D = = (x 2 +7n - (x 1 + 712) — ) yi + (x 4 + 713)2/3 + (x 3 + 714)2/4 

v Pl2 y 

+(x\ + 7i 2 )/3n + 711^1 + 7i2^ 2 + 713X3 + 714X4 + 711712 + 713714 • (5.10) 
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After getting three cut equations of four variables, we can use the Grobner basis method with ordering 
(x4,yi,y4,y3) and constraints on the powers of monomial x^ X4 ^ 1 y^ V3 ^ 'yf^ 

E d( Xi )<5, ^)< 3 ' E E d (^ 6 ( 5 - n ) 

all ISPs of x all ISPs of y all ISPs of x all ISPs of y 

to find integrand basis. Elements of integrand basis will be different depending on actual kinematic 
configurations. In this case, there are three kinds of integrand basis depending on whether if K4 is massless 
or K$,Kq is absent. For all kinematic configurations with K4 massive, the integrand basis contains 14 
elements given by 

#542 = { 1 ,x4,yi,y3,x4y3,yiy3,yhyhy4:,y3y4,y'iy4:,yl,y3yl,yl} ■ (5.12) 

For kinematic configurations with if 4 massless while at most one of K^,Kq to be absent, the integrand 
basis contains 14 elements given by 

Bb42 = i l iXiiyi,y3,X4yz,y 1 yz,yl,y 1 yl,yl,y4,y 1 y i ,yl,y 1 yl,yl} . (5.13) 

For kinematic configurations with K4 massless and both K§ = K§ = 0, the integrand basis has 20 elements 
given by 

#B42 = { 1 ,x4,yi,x4yi,yl,x4yl,yl,x 4 yf,y3,yiy3, 

yfy3, yl,ywl,yi, V4,yw4, yfy4, yl,ywl, yl) ■ (5.14) 

Note that elements of integrand basis generated from different ordering of ISPs will possibly be different, 
but after choosing one ordering, there will always be three kinds of integrand basis depending on kinematic 
configurations. 

After given integrand basis, we need to discuss how to get their coefficients from integrand input 
generating from Feynman diagrams or unitarity cut method. As in the one-loop case, either algebraic 
geometry method or parametrization method can be used. 

The algebraic geometry method is illustrated as follows. Firstly we should get integrand expression 
input T{l\,l2) from Feynman diagrams or unitarity cut method after subtracting contributions from 
higher topologies. After expanding £1,^2 into momentum basis and substituting all kinematics with ISPs, 
we can rewrite T{l\,l<i) as polynomials of ISPs, in this example J r (x4, y±, 2/3, 2/4). Then we can divide 
^(x4,yi,ys,y4) by Grobner basis generated from ideal / = (^Dq, Do, Dq^ with ordering (X4, y±, 2/4, 3/3). 
The remainder of division is linear combinations of all terms in integrand basis with wanted coefficients. 

All coefficients can be found at the same time by using above algebraic geometry method, but it 
may take long time to do so if the number of elements is large. Instead we can use branch-by-branch 
polynomial fitting method (see reference [29]) to simplify problem, and fit a smaller set of coefficients 
at one time. The idea can be illustrated as follows. Because Dq = a(x 4 — x^ 1 )(x 4 — x\ 2 ), we can divide 
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polynomials F(xi, y±, y 3 , y 4 ) by Grobner basis generated from I\ = <^(x 4 — x A 1 ),Do,DoJ with ISPs ordering 
(#4, yi, y4, y 3 ). After the division T{x±,y\,y$,y4)/I\, we will get remainder 

/ R{J 7 (x4,y 1 ,y 3 ,y 4 )/I 1 ) = fi + / 2 y3 + hvl + Uvl + hVA + hvzVA + hvh& , ( 5 - 15 ) 

with seven known coefficients /j. It is easy to see that the remainder of 14 integrand basis over I\ is given 
by 

{l)/h ->• 1 , (x 4 )/il ->• d 2 , {yi)/h ->• ^31^4 + fey3 + fe , (?/3)/A -> ^3 , 

(xm)/h ->• d 5?/3 , (ym)/h ->• fey3y4 + fey! + feys , ->• y! 5 (yf)/ 7 ! ->• yf » 
(y4)/A ->• y4 , {yzyi)/h ->• y3y4 , (yhd/h -> yfy4 , 
(yD/A ->• du^myA + fe,2y4 + di2,3yi + rf i2,4y3 + fe,5 , 
{yzyl)/h ->• ^i3,iy!y4 + di 3i2 y 3 y4 + d^yi + ^i3,4y! + di3,sy3 , 

(yD/A di4,iy!y4 + di^ 2 yzy± + di4,3y4 + du^yf + di4,sy! + ^i4,6y3 + di4,7 , (5.16) 

with known coefficients c2. Thus by comparing both sides we obtain following seven equations of 14 unknown 
coefficients Cj from one branch 

fl = Cl + C2(i2 + c3fe + Cl2fe,5 + cuduj , 

h = c 3fe + C 4 + C 5 d 5 + C 6 fe + Cl 2 dl2,4 + Ci 3 di 3j 5 + Ci 4 di4 i6 , 

/3 = Qjfe + C 7 + Ci 2 di2,3 + Cl 3 di 3i 4 + Ci 4 di 4j5 , 

/ 4 = C$ + 013^13,3 + Ci 4 di 4j4 , 

h = c 3 d 51 +Cg+ Ci 2 dl2,2 + Ci 4 di 4j3 , 

h = ced 6 i + cio + ci 2 fe,i + ci 3 di 3i 2 + cudu :2 , 

/7 = en + ci 3 di 3) i + c M d Ut i . (5.17) 

Similarly, we can divide polynomials -F(x 4 , y±, y 3 , y 4 ) by Grobner basis generated from another branch 
l 2 = ((x& — X4 2 ), Dq,DqJ with the same ISPs ordering (x 4 , y 4 , y 4 , y 3 ). After that we can get another seven 

equations relating to q with other known coefficients d. With this modified algebraic method, we can 
get a smaller set of coefficients in one branch, e.g., in this example each branch can be used to write down 
seven equations (we will say that this branch can detect seven coefficients). Combining results of both 
branches we get 14 independent equations, and they can be used to solve 14 coefficients of q. 

Besides algebraic geometry method, it is also possible to find coefficients by parametrization method. 
This method is tightly related to the branch-by-branch fitting method. In this example, we can use -Do to 
solve X4 and get two solutions. Then we put one solution x 4 ' to Do, Do, and use one variable, for example, 
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U4 to express 2/1,2/3- Finally we put 2/1(2/4), 2/3(2/4) back to the identity 



14 

F{x±\ 2/1(2/4), 2/3(2/4), 2/4) = ^c k B B A2,k (2/4) , (5.18) 

fc=i 

and find coefficients Cj by comparing both sides. This method is very useful to evaluate coefficients 
analytically or numerically. In this example, we only need to take arbitrary 7 values of 2/4 to produce 
seven equations from each branch, and solve 14 linear equations by combining two branches to find all 
coefficients. 

5.1 Structure of variety under various kinematic configurations 

In some kinematic regions, for instance, some of external momenta being massless or absent, the variety 

will split into different branches. Analysis of branch structures is helpful to find coefficients of integrand 

basis in actual calculation. Notice that in this example there are four equations containing only £1, and 

r r 

by solving these equations we can fix X4 to two solutions x 4 1 ,x 4 2 , no matter what the external kinematic 
configuration is. Thus we only need to focus on the two remaining equations Do, Do with X4 replaced 
by two solutions x 4 1 ,x 4 2 . Since in general x 4 * / x 4 2 , branches parameterized by different x 4 * will not 
intersect with each other. 

When K4 is massive, /?n 7^ 0, the on-shell equation Do = is non-degenerate. If we take 2/1 = t as 
free parameter, D = becomes non-degenerate conic section of variables 2/3,2/4 while D = becomes 
linear equation of variables 2/3,2/4- Using following two equations 

D = : y 3 y 4 + F(r) = , A) = : ay 3 + by 4 + c(r) = , (5.19) 

where a, b are some constants, F(t) is second order function of r, and c(r) is linear function of r, we can 
solve 

ac{T)±^a^c{Tf + AabF{r)} 
2/4 = . (5.20) 

2/4 is a rational function of r if c(r) 2 + AabF(r) inside the square root is a perfect square of r. Using the 
explicit expression of F(t), c(t) and a, b we find the discriminant of quadratic function c(r) 2 + 4a6F(r) to 
be 

{x\xl + xlxlWn (4 + 7ii + PuM + M + ^ + (*4 + TiaM + 714) g 

P12 P12 

where 

2 = 711X1 + 712X2 + 713X3 + 7ux 4 + 711712 + 713714 , (5.22) 

and x\ denotes the solution of x\ with X4 = x 4 . The first term in above result vanishes because Dq = 
%i x 2 + X3XJ = 0. Generally the second term will not be zero, but if = 0, i.e., 7h = 0, we have H = 0. 
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Similarly, if K 5 ^ but K & = 0, using 711 = -j3 n - a 31 , 712 = -P12 - «32, 7i3 = -«33 and 714 = -a 34 , 
the second term becomes 

-x\azi - x\ a 32 - x^a 33 - x\au + a 3 ia 32 + 033034 n _ £>3 - A) CT ^ 



£12 2/3i 2 

which vanishes because of the on-shell equation D 3 — Dq = 0. In these cases, c(r) 2 + AabF(r) is a perfect 
square, and we can get two solutions which are rational functions of free parameter y\ for each solution 
x\. In other words, each irreducible branch will split into two branches in these special kinematic region, 
and in total we get four branches denoted by V^ l > l , y r i> n 2 anc | yr 2 ,rii^ yr 2 ,n 2 _ Yiaich. branch can detect 
4 coefficients, and the two branches V^ 1 ^ 1 ,V Tl ^ 2 intersect at a single point. Similarly, the two branches 
V r ' 2,Tl1 , y r2 > n 2 intersect at another point. There are no intersection among other combination of branches. 

If is massless, i.e., /?n = 0, but at most one of K§, Kq is absent, then Dq = = 2/32/4. There are two 
branches parameterized by 2/3 = with 2/4 free parameter or 7/4 = with 2/3 free parameter. Considering the 
remaining linear equation Dq = of (2/1,2/3, 2/4), it is easy to see there are also four branches V ri,ni , V Vl ' U2 
and y 1 " 2 ' 111 , y r 2> n 2 The intersection pattern of these four branches is the same as in previous paragraph 8 . 

If we further consider a special kinematic region where K4 is massless and both K$,Kq are absent, 
unlike previous discussions, the dimension of variety will increase from one to two, and the integrand basis 
is given by 20 elements as shown in ( 5. J 4 ) instead of 14 terms. This can be explained by noticing that the 
three equations become 

«24 2 , «23«24 ~ ai2«21 + "21022 

A) = x A H x 4 , Dq = y 3 y 4 , 

«23 «23 

f. "24 . a21«22 - "21012 + «23«24 n . s 

Do = ^42/3 ^42/4 H 2/4 (5.24) 

"23 a 2 3 

in this special kinematic region, and 2/1 disappears. The variety is given by two branches. One branch is 
parameterized by X4 = 0,2/4 = while 2/1 , 2/3 f ree parameters, and the other branch is parameterized by 
X4 = (0210:22 — CK12CK21 + a23«24)/a24 5 2/3 = while 2/1,2/4 free parameters. The degree of freedoms has 
increased from one to two. Each branch can detect 10 coefficients and there is no intersection between 
them, so adding them up we can detect all 20 coefficients. 

6. Example two: non-planar crossed double-triangle 

Our second example will be non-planar crossed double-triangle topology (C22) as shown in Figure |6|. Differ- 
ent from planar penta-triangle topology (B42), the variety of (C22) is two-dimensional, so the intersection 
between different branches would be one-dimensional variety instead of single points. Topology (C22) also 
has symmetry of relabeling (Ki, K 2 , K 3 ) as well as symmetry of relabeling (if 4, K$). Discussion of different 
kinematic configurations can be simplified by using these symmetries. 



Besides branch structure of variety, the integrand basis (5.12) need to be modified too. The reason is that in the case 



Kl = 0, we have Do = 3/33/4, thus elements such as J/3 J/4, J/3 J/4, J/3 J/4 could be divided by Dq and they are no longer elements 



of integrand basis. The modified integrand basis is given by (5. IS). 
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(C42) ( C33) (C32) (C22) 

Figure 6: The type (C) contains all 4 diagrams with n 3 = 2. Every diagram is denoted by (Cnm) where n,m are 
the numbers of propagators containing pure l\ or £2 loop momentum. The diagrams are drawn in most general form, 
and some external momenta, for instance Kq,Ky in (C42), could be absent. All external momenta are out-going 
while convention of loop momenta is labeled as an arrow in each diagram. 

For this topology, we use K±,K 2 to construct momentum basis and expand external and loop 
momenta as 

4 4 
K\ = anei + a 12 e 2 , K 2 = fi\\e\ + /3i 2 e 2 , if 4 = ^ 71^ , K 3 + K A = ^ 7 2i e; , 

i=l i=l 

h = x 2 ei + xie 2 + x 4 e 3 + x 3 e 4 , t 2 = y 2 e x + y±e 2 + y 4 e 3 + 7/ 3 e 4 . (6.1) 

With this expansion, 6 propagators can be rewritten as functions of 8 RSPs. The two propagators that 
containing only x% variables are given by 

Do = tf = 2(x\x 2 + x 3 x 4 ) , Di = (£1 - Ki) 2 = D + 2a\\a\ 2 - 2{a\\X\ + ai 2 x 2 ) . (6.2) 

The two propagators that containing only j/j variables are given by 

D = £\ = 2( ym + ym ) , Di = (£ 2 - K 2 f = Do + 2/3 n /3 12 - 2(p im + /3 12 y 2 ) . (6.3) 

The remaining two propagators contain both Xi,yi variables 

Do = {h + i 2 + K 4 ) 2 = D0 + D0 + 2(7n7i 2 + 713714) + 2(xiy 2 + x 2 y x + x 3 y 4 + x 4 y 3 ) 

+ 2(7n(xi + yi) + 7i 2 (x 2 + y 2 ) + 7i 3 (x 3 + y 3 ) + 7i 4 (x 4 + y 4 )) , 
D\ = (h + £ 2 + K 4 + K 3 ) 2 = D + D + 2(7 2i 7 22 + 7 23 7 24 ) + 2{x x y 2 + x 2 yi + x 3 y 4 + x 4 y 3 ) 

+ 2(7 2i (xi + yi) + 7 22 (x 2 + y 2 ) + 7 23 (x 3 + y 3 ) + 7 24 (x 4 + y 4 )) . (6.4) 

Three linear equations D\ — Do = 0, D\ — Do = and D\ — Dq = can be constructed from 6 on-shell cut 
equations, and from which we can solve xi,y 2 ,x 2 as functions of five ISPs (x 3 , x 4 , y±, y 3 , y 4 ). Substituting 
these solutions back into Do, Do, Do we get three polynomial equations, which define the variety of this 
topology. 

We will consider the kinematic configuration where K4, K§ is absent, or some of K\, K 2 , K 3 are mass- 
less. In order to make the kinematic configuration clear, we use the notation C22^fjpl R \ where each L, N, R 
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could be either M or m representing massive or massless limit of K\,K 3 ,K2 respectively. U, P could be 
either K A ,K§ if they are non-zero or if they are absent. In this notation, for example, 022^^^ 
represents kinematic configuration with K\,K 3 massive, Ki massless, K A non-zero and K5 absent. 



6.1 The integrand basis 

To deti 
ditions 



To determine the integrand basis, we take all possible monomials x^^yf^ y*^ under con- 



d( Xi )<4, d (^)<4, E d (^)<5, (6.5) 

all ISPs of x all ISPs of y all ISPs of x all ISPs of y 

and divide them by Grobner basis generated from Do, Do, Do with ISPs ordering X3 > y 3 > x A > y A > y\ 
(we will use the same ordering through this example) and collect all monomials in the remainder. For 
different kinematic configurations, the number and elements of integrand basis can be different. The 
reason for this can be understood by checking the example ay\ + 1/37/4. If a ^= 0, 2/32/4 can be taken as 
integrand basis, but if a = 0, J/3J/4 is reducible and can never be integrand basis. 

After checking all 24 different kinematic configurations, we find that there are in total 6 different kinds 
of integrand basis. For kinematic configurations where at least one of K A ,K§ are non-zero and K\,Ki are 
massive, the integrand basis contains 100 elements given by 

{1, x 3 , x A , x 3 x A , x\, X3X4, x% x 3 x\, x\, y\,x 3 y\,x A y\,x 3 x A y\, x\y\ , x 3 x\y\ , x\y x , x>,x\y\ , x^yi , yf, 
x 3 y\,x A y\,x\y\, x\y\,y\,x 3 y\,x A y\,x\y\,y\, x 3 yf, x A y\, y\, y 3 , x 3 y 3 , x A y 3 , x\y 3 , x\y 3 , x\y 3 , yiy 3 , 
x^yiys, x\yiy 3 , x\yiy 3 , y\y 3 , x A y\y 3 ,y\y 3 , x A y\y 3 , y\y 3 ,yl,x 3 yl, x A y 3 ,yiy 3 ,x A yiy 3 ,yly 3 ,x A yfy 3 , 
ylyhyl,x 3 yl, x A yl,y x yl,x A yiyl, y\yl, yf, x 3 y|, x A y 3 , y x y\, y\,y A ,x A y A , x\y A , x\y A , x\y A , y x y A , 
x A yiy A , x\y x y A , x\yyy A , y\y A , x A y\y A , x\y\y A , y\y A , x A y\y A , y\y A , y\,x A y\,x\y\, x\y\,y x y\,x A y\y\, 

2 2 2 2 22323 323 3 3234 4 4 5i I a a\ 

x A yiy A ,y 1 y A ,x A y 1 y A ,y 1 y A ,y A ,x A y A ,x A y A ,yiy A ,x A y 1 y A ,y 1 y A ,y A ,x A y A ,yiy A ,y A \ . (6.6) 

For kinematic configurations where at least one of K A ,K§ are non-zero and K\ is massive, K% is massless, 
the integrand basis still contains 100 elements, and is given by replacing one element from ( |6.6D 

#C22 = #C22 " {4vl} + {xWm} ■ (6.7) 

For kinematic configurations where at least one of K A ,K^ are non-zero and K\ is massless, the integrand 
basis contains 98 elements, and is given by removing 17 elements from (|6.6| ) while adding another 15 
elements: 

#C22 = ^c-22 - {^32:4, x 3 x\,x 3 x\, x 3 x A yi,x 3 x\yi,x 3 x\yi,x 2 A y 3 , x\y 3 ,x\y 3 ,x\yiy 3 ,x\yiy 3 ,yl, 

2 3 22 2 2 5i , r 2 3 4 2 3 4 2232232 

x A yiy A , x A yiy A , x A y x y A , x A y x y A ,y A \ + {x 3 ,x 3 ,x 3 ,x 3 y 1 ,x 3 yi,x 3 y 1 , x^x^x^x^, 
z|3/3, xjy 3 , xlyl,xlyl,xlyl) . (6.8) 
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If both K^,K§ are absent and K2 is massive, the integrand basis contains 96 elements, and is given by 
removing 22 elements from (|6.6| ) while adding another 18 elements 



&C22 = &C22 ~ {^3x4, x 3 xj, x 3 x% x 3 x i y 1 ,X3xly 1 ,x 3 xly 1 ,x 3 yj, x\y{,x\y\, x 3 yf,xjyf, x 3 yf, 
xly3,xly 3 ,xiy 3 ,xlyiy 3 ,xlyiy 3 ,yiyi, y\,x\y\y^ yxyj, y\} + {x\,x\,x\, x\yx,x\yx, 
xiyx,x 2 3 y 3 ,xly 3 ,xly 3 ,x 3 yxy 3 ,x 2 3 yiy 3 , x\yxy 3 , x\y\, x\y\ , x 3 yxyl,xlyxyl,xlyl,x 3 yxyl] . (6.9) 

If both K^,K^ are absent, if 2 is massless, and at least one of Kx,K 3 are massive, the integrand basis 
contains 96 elements, and is given by replacing 9 elements from (|6.9| ) 



nV nlV ( 2 2 2 2 3 342 3 2 2i 

, r 2 22 32 22 32 32323 4t 1 n\ 

+{^32/i, x 3 y x , x 3 2/i, x 4 y 1 ,x 4 y 1 , x 3 y Xl x 3 y 1 ,x 4 y 1 , x 3 y x \ . (6.10) 

Finally if both K^K^ are absent, and all Kx,K 2 ,K 3 are massless, the integrand basis contains 144 terms, 
which is given by 

aVI 
°C22 — 

r, 234 232222334 234 

\x, xx, Xi, Xi, Xi, x 2 , xxX2, x x X2, X1X2, x 2 , xxx<ii X1X2, X2, xxx 2 , x 2 , 2/1, xxyxt x x yx, ^i2/i, x x yx, 
X2y\,xxX2yx,x\x2yx,x\x 2 yx,x\yx,xxx\yx,x\x\yx,x\yx,xxx\yx,x\yx,y\, xxy\,x\y\, x\y\, 
X2vl,xix 2 yl, x\x 2 y\, x\y\,xxx\y\, x\y\,y\, x x y\, xjyf,x 2 yf, x x x 2 y\, xlyf,yf,xxyf, x 2 yf, 
2/2, xxy2, xjy 2 , xfy 2 , xfy 2 , x 2 2/2, xxx 2 y 2 ,x\x 2 y 2 , x\x 2 y 2 ,x\y 2 ,xxx\y 2 , x\x\y2,x\y2,xxx\y2, 
X2V2, 2/12/2, xxyxy2, xjyxV2, xfyxV2, ^22/12/2, ^1^22/12/2, xjx 2 yxy 2 , x\yxy 2 , ^1^2/12/2, £22/12/2, 
2/12/2, xxy\y 2 yx\y\y 2 , x 2 yfy 2 , x x x 2 y\y 2 , x\y\y 2y y\y 2 , xxyfy 2 , x 2 yfy 2 , y\, xxy\,x\y\, x\y\, 
X2vl,xxx 2 yl, x\x2y\, x\y\,xxx\y\,x\y\,yxy\,xxyxy\, x\yxy\,x 2 yxy\-,xxx 2 yxy\, x\yxy\-, 

2/12/2, £12/12/2, £22/12/2, yh x iyh x iyt x 2yh x i x 2yt x 2vt 2/12/I, £12/12/2, £22/12/2, 2/2, £12/2, £22/2, 

2/3, ^12/3, £i2/3, £?2/3, £i2/3, 2/12/3, £l2/l2/3, £i2/l2/3, £?2/l2/3, 2/i2/3, £i2/i2/3, £i2/i2/3, 2/i2/3, £l2/i2/3, 2/i, 

22232 2 22 222 223 323 3 34 4t /cn\ 

xiys^iy^x^yxysjXiyxya^^xy^y^xxyiysjys^xy^x^yiy^xxyiys^sjXiys} . (6.11) 
6.2 Structure of variety under various kinematic configurations 

Having given the integrand basis we move to the discussion of variety determined by six propagators under 
various kinematic configurations. 

6.2.1 Kinematic configurations with K 4 ,K^ non-zero 

Given the integrand basis, the problem then becomes fitting coefficients of integrand basis for given diagram. 
As mentioned above, the computation can be simplified when using branch-by-branch method, thus it is 
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important to study the structure of variety in various kinematic configurations. For general case where both 
K4,Ks non-zero and Ki, K2, K% massive, the variety defined by six on-shell cut equations is irreducible, 
i.e., there is only one branch with dimension two. All 100 coefficients of integrand basis ( |6.6| ) should be 
determined at the same time using this irreducible branch. 

The variety will split into two branches when one of K\, K2, K% is massless, this corresponds to kine- 
matic configurations C22^^™\ 022^"^^ ■* and C22^'^^\ It is easy to see that when K\ = 0, fin = 
so D = ysUi, or when Kf = 0, a\2 = and D = X3X4. For = 0, We could use the massless condition 
(721 - 711) (722 - 712) + (723 - 713X724 - 7u) = to solve 724, and substitute it back to D ,D to solve 
2/3, £4. After putting solutions of 2/3,2:3 back to Dq, the numerator of -Do is factorized into two factors, 
which gives two branches. Above procedure, although can be applied to both numeric and analytic modes, 
is straightforward but complicated, and probably miss some branches in certain kinematic configurations. 
An alternative and better way of finding branches of variety is to use Macaulay2 [36] , which can systemat- 
ically get the information of branches. Let us take kinematic configuration C22^^^ as an example to 
illustrate the structure of these two branches. In this example, one branch is characterized by 2/3 = and 
the other branch by 2/4 = 0. For the first branch 2/3 = 0, only 65 elements are left after assuming 2/3 = 
for all elements in ( |6.7D . Dividing these 65 monomials over Grobner basis generated from cut equations 
defining this branch, we find that only 59 of them are independent. So we can fit only 59 coefficients of 
integrand basis fl6,7| ). Similarly, for the second branch 2/4 = 0, 66 elements are left after assuming 2/4 = 0, 
and only 59 are independent after dividing 66 monomials by Grobner basis generated from three cut equa- 
tions of this branch. Both branches are varieties of dimension two, and their intersection is an irreducible 
variety of dimension one. By using the intersection one can detect 18 coefficients, which means that we 
can find all 59 + 59 — 18 = 100 coefficients when using both branches. 

If two of K\, K 2 , K$ are massless, i.e., kinematic configurations C22^"^ , C22^'^ n ^ > and C22|™^™\ 

the variety will further split into 4 branches. Take kinematic configuration C22^^.^ as an example, mass- 
less conditions of K\, K2 will reduce Dq = X3X4 and D$ = 2/32/4- It is easy to see that there are 4 branches 
characterized by V\ : (x^ = 0, 2/3 = 0), Vi : (a?3 = 0, 2/4 = 0), V3 : (X4 = 0, 2/3 = 0) and V4 : {x± = 0, 2/4 = 0). 
Using algebraic or other methods, one can find that each branch can detect 34 coefficients. A naive summa- 
tion of these 4 branches gives 34 x 4 = 136 coefficients, which is larger than the number of integrand basis. 
This means that there are intersections among 4 branches. By analyzing intersections among all possible 
combinations of branches, we find that intersections for pairs (Vi, V2), (V±, V3), (V4, V3), (V4, V2) are irre- 
ducible one-dimensional varieties 9 , and for each pair (^1,14) and (V2,Vs) isolated points. For intersection 
among three or four branches, it is located at above two isolated points. 

If we assume kinematic configuration to be C22|^"^™y ) where all K\, K2, K3 are massless, the variety 



9 Each one-dimensional intersection can detect 10 coefficients for this example. With information of other intersections, 
we can make following counting. Since intersection of three or four branches detects 2 coefficients, each intersection of two 
branches will detect 10 — 2 = 8 independent coefficients, thus each branch will independently detect 34 — 8 — 8 — 2 = 16 
coefficients that can not be detected by other branches. Adding all together we have 16x4 + 8x4 + 2 = 98 coefficients as it 
should be. 
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is given by eight branches, i.e., each branch of previous paragraph has further split into two branches. 
The first two branches Vi,V 2 characterized by X3 = 2/3 = (or the seventh and eighth branches Vj,V& 
characterized by x 4 = 7/4 = 0) can detect 19 and 21 coefficients respectively, and 34 coefficients can 
be detected by using two branches. This can be checked by noticing that the intersection of these two 
branches can detect 6 coefficients, so 19 + 21 — 6 = 34. Similarly, each of the third and fourth branches 
V3, V4 characterized by X3 = j/4 = (or the fifth and sixth branches V5, Vq characterized by x 4 = 2/3 = 0) 
can detect 20 coefficients, and 34 coefficients can be detected by using two branches. This can also 
be checked by noticing that their intersection can detect 6 coefficients, so 20 + 20 — 6 = 34. We also 
need to clarify the intersection pattern among eight branches. There are no intersections shared by five 
or more branches. The intersection of following six pairs (Vi, V2, V3, V4), (V5, Vq, V?, V$), (V±, V 2 , V5, Vq), 
(V3, V4, V7, Vs), (Vi, V3, Vq, V$), (V2, V4, V5, V7) is single point. Intersection of every three branches is also 
single point, which is inherited from corresponding intersection of every four branches, such as intersection 
point of (Vi,V2,V"3) coming from intersection point of (Vi, V2, V3, V4). No new intersecting point besides 
the ones of every four branches is found in intersections of every three branches. The intersections of every 
two branches are possibly one-dimensional varieties or single points. In order to express the intersection 
pattern, we will use following notation V\ fl V2 = (d\m) where d is the dimension of variety (so d = 1 for 
one-dimension and d = for points) and m is the number of coefficients detected by the intersection. Thus 
all possible intersections between pair (Vi, Vj) are given by 



(1|6) 


= v 1 n v 2 


= v 2 n v 4 


= V 2 n Vq 


= v 3 n v 4 


= v 3 n v 8 


= V 5 n Vq 


= v 6 n v 8 


= v 7 n v 8 


(1|5) 


= Vi n v 3 


= v 1 nv 6 


= v 4 n v 7 


= v 5 n v 7 










(0|1) 


= Vi n v 4 


= v 1 nv 5 


= Vx n v 8 


= v 2 n v 3 


= v 2 n v 7 


= V 3 n Vq 


= v 3 n v 7 


= v 4 n v 5 




= v 4 n v 8 


= v 5 n v 8 


= Vq n v 7 













6.2.2 Kinematic configurations with one of K±,K§ absent 

For C22x 4) or C22k 5 ,0, i.e., one of K 4 ,K^ absent, the variety is given by two branches 10 even without 
imposing massless conditions of Ki,i = 1,2,3. Each branch can detect 64 coefficients of integrand basis 
and their one-dimension intersection can detect 28 coefficients. By using two branches all 64+64 — 28 = 100 
coefficients can be detected. 

For kinematic configurations C22™^, M , C22^^^ and C22^^™ , the variety is given by four 
branches. To illustrate the structure of each branches, let us take C22^5^ m as an example. Each branch is 
2-dimensional variety and can detect 21 coefficients. Let us use Vi, V 2 to denote two branches characterized 
by V3 = 0, and V3, V 4 to denote two branches characterized by y 4 = 0. We find that these four branches 
will intersect at a single point. Among intersections of every three branches, non-trivial two intersecting 
points exist for pair (Vi,V 2 ,V 4 ) and (V 2 , V3,V 4 ). Pair (Vl,V3) intersects at a point, while intersections of 

10 This can be seen by solving 1/3, X3 using Do = and D = equations and putting solutions back to Do, which is factorized 
to two pieces. One can also use Macaulay2 to find branches. From now on, we will not discuss how to get branches. 
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all other five pairs of every two branches are one-dimension varieties. Among them V± D V2, V3 fl V4 can 
detect 11 coefficients while V\ n V4, V2 n V3 can detect 6 coefficients and V2 n V4 10 coefficients. It is also 
worth to mention that though having four branches, the intersection pattern of C22™^^ , C22^^^ 
and C22£$£ are different from C22™;>f,C22* and C22^- 

Next let us discuss kinematic configurations C22^™^ , C22^^T™ and C22^^^ . For these 
cases, the variety is given by six branches. Taking C22^^™ as an example, the first two branches V\, V2 
characterized by X3 = ys = can detect 19 and 21 coefficients respectively, and the intersection of these two 
branches can detect 6 coefficients, thus we have 34 coefficients by using both branches. The third branch 
V3 characterized by Z3 = j/4 = can detect 34 coefficients. Similarly, the fourth branch V4 characterized 
by X4 = 2/3 = can also detect 34 coefficients. The last two branches V5, V% characterized by X4 = 1/4 = 
can detect 19 and 21 coefficients respectively, and by using both of them one can detect 34 coefficients. 
It is interesting to notice that these six branches are split from corresponding 4 branches of 022^^™. 
We will again clarify the intersection pattern of these six branches. No intersections exist for every five 
or six branches. For intersection of every four branches, pair (Vi, V3, V4, V5) and (V2, V3, V4, Vq) intersect 
at single points. Apart from the inherit intersecting points of four branches, there are also pairs of every 
three branches (Vi, V2, V3), (V4, Vj>, Vq), (Vi, V2, V4), (V3, V5, Vq) that intersect at different single points. For 
intersection of every two branches, (V\, V5), (V2,Vq) intersect at one single points, and (V3, V4) intersects 
at two points, and (V-j^V^), (^l,^), (V§,Vq) intersect at one-dimensional variety which can detect 6 
coefficients, while (^3,4, ^1,5) also intersect at one-dimensional variety which can detect 5 coefficients. 

For kinematic configuration C22™^J™ where all Ki,K2,Ks are massless, the variety splits to eight 
branches. The branch structure is the same as C22^"^-™. Two branches V\,V2 characterized by Z3 = 
2/3 = as well as two branches V7, Vs characterized by X4 = 1/4 = can detect 19 and 21 coefficients 
respectively, while two branches V3, V4 characterized by X3 = 1/4 = and two branches V5, V% characterized 
by X4 = 7/3 = can detect 20 coefficients respectively. These eight branches intersect at a single point, while 
all intersections among every seven, six, five, four or three branches are also located at the same point. 
There are 28 possible pairs of intersecting two branches, among them 12 are one-dimensional varieties, 
and intersections of the remaining 16 pairs are the same single point as intersection of eight branches. 
For the one-dimensional variety, 8 of which coming from (^2,^1,4,5), (^,^6,7), (V3,V<0 and (V$,Vq) can 
detect 6 coefficients individually, while the other four that coming from (Vi, Vjj^) and (V7, V^s) can detect 
5 coefficients. 



6.2.3 Kinematic configurations with both K4,K$ absent 

For kinematic configuration C22 0) 0, since K4 = = 0, momentum conservation ensures K3 = —K\ — K2. 
K±,K2 are still independent, so we can use them to construct momentum basis e^. For this simple case, 
we can write down analytic expressions and make discussion more transparent. 

Using parametrization K\ = a.\\e\ + a\2&2 and K2 = + Pvie-i-, the three non-linear cut equations 
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can be given by 

011012 yi > 
D = X3X4 H — 3 — (l - -3— )yi , 
P12 P12 

-Do = 2/32/4 + Pn(l - -5 — )2/i , 
Pi 2 

-Co = X4J/3 + H (1 - — )y 1 (6.12) 

P12 P12 

after eliminating all RSPs. If K\,K2,K 3 are massive, the variety is given by following six branches defined 
by ideals: 

™„(M,M,M) ™„(M,M,M) 

-°^(0,0) _ r„, „, \ \r U (0,0) 



V l (0 - 0) = {2/3,^,2/1}, V 2 ^ ={y 3 ,x 3 , yi -f3 12 } , 

T7 • (0,0) r t t 7" (0,0) r a t 

^3 = {2/4,2/1,^4} , V4 = {2/4,2/1 - £12,3:4} , 

C22 (M,M,M) 

V 5 (0 ' 0) = {y 3 2/4 + £n(l - 2/i/£i2)2/i,2/3ai2 - x 3 £i2, 2/4aii - aW^ii} , 

C22 (M,M,M) 

^ 6 (0 ' 0) = {2/32/4 + £ii(l - 2/i/£i2)2/i, -2/3aii + a:3£ii,2/4Q!i2 - x 4 /3i 2 } ■ (6.13) 

Among these six branches, four of them = 1,2,3,4 will individually detect 19 coefficients and two of 
them VI, i = 5, 6, 36 coefficients. The physical picture is following. Each branch of C22^^j^ will split 
into three branches with two branches detecting 19 coefficients and one branch detecting 36 coefficients. 
The intersection pattern of six branches are following. No intersection exists for six or every five branches. 
Each combination of (V2, V4, V5, Vq) and (Vi, V3, V5, Vq) intersects at a single point. All intersection of 
every three branches are coinciding with the single point of four branche intersection, and there are no new 
intersections among other combinations. For intersections of every two branches, there is no intersection 
among four pairs (^1,3,^,4), and (Vi,V 3 ) intersects at one point, while (V^VJj) intersects at two points. 
Intersections of remaining nine pairs are one-dimensional variety. 

If one or two momenta of Ki, K2,K 3 are massless, i.e., kinematic configurations C22^^ M \ C22^ "' W 



'(0,0) ' (0,0) 

C22^'^' m \ C22^'™j M \ C22|™'^'"^ and C22^'^' n \ the variety still has six branches. Definition of 



branches are still the same as given in ( S.13| ) for C22^™'^i' M \ C22^'™? M \ C22^'™^ M \ but will be dif- 



(0,0) ' (0,0) ' (0,0) 



ferent for C22^'^' m \ C22|™'*^' m \ C22^'^' m ^ where the first four branches keeps the same as given in 



( |6.13 ), and the last two branches are given by the ideals 



Q^t^M ,M ,m) ,{m,M ,m) ,(M ,m : m) 

V 5 (0 ' 0) = {?/3, -2/4«12 + X A (5i2,X 3yi + Qn (1 — 2/l/>5i2)2/l} , 

V 6 (0 ' 0) = {2/4,2:42/3 + an (l-2/i /£i2)2/i, -2/3"i2 + X3012} . (6.14) 

For the last kinematic configuration C22fc^t' m ^ , the external momenta are extremely degenerate since 
we must either have Ai ~ A2 ~ A3 or Ai ~ A2 ~ A3. In other words, we can not use K%,K2 to construct 
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momentum basis ej. One possible choice of momentum basis is the massless momenta K\,K 2 , e 3 ,e 4 
satisfying 11 

Ki ■ K 2 = K\ ■ e 3 = K 2 ■ e 4 = e 3 • e 4 = , K\ • e 4 = K 2 ■ e 3 = 1 . (6.15) 

With this momentum basis we can expand loop momentum t\ as 

h = {h • e A )Ki + (h • e 3 )K 2 + {h ■ K 2 )e 3 + {h ■ Ki)e A = xiKx + x 2 K 2 + x 3 e 3 + x 4 e 4 , (6.16) 

and similarly for l 2 . Then the six propagators are given by 

D = £j = 2{ Xl x 4 + x 2 x 3 ) , D 1 = (h - K\) 2 = D - 2x 4 , 
D = l\ = 2( yi y 4 + y 2 y 3 ) , B x = (£ 2 - K 2 f = D - 2y 4 , 
D = {h + h) 2 = 2(xi + yi )(x 4 + Va) + 2(x 2 + y 2 )(x 3 + y 3 ) , 

Di = (ii + h + K 3 ) 2 = D - 2(x 4 + y 4 ) - 2(x 3 + y 3 ) . (6.17) 

We can easily solve these equations and get 

x 4 = , y 4 = , x 3 = -y 3 , (6.18) 

and there are only two non-linear equations left 

D = x 2 y 3 , D = y 2 y 3 . (6.19) 

Now we have five ISPs (xi,x 2 ,yi,y 2 ,y 3 ) and two non-linear equations expressed with these ISPs. The 
integrand basis can be given by monomials of ISPs, and there are 144 elements under degree conditions as 
already shown. The variety is given by two branches. The first one characterized by y 3 = is dimension four 
variety, and the second one characterized by x 2 = y 2 = is dimension three. The first branch can detect 
114 coefficients while the second branch can detect 49 coefficients. Their intersection is two-dimensional 
variety and detects 19 coefficients. 

The splitting of branches in different kinematic configurations is summarized in Figure [?]. 

7. Remaining topologies 

After demonstrating methods and various properties with above planar penta-triangle and non-planar 
crossed double-triangle examples, we will present results for the remaining topologies in this section. We 
will omit many details but just show the main results for each topology since analysis is quite the same as 
previous discussions. 

n We can always have this choice. For example, if K\ = AiAi, K2 = A1A2 we can take e% — C3A2A1 and en = C4A2A2. 
Similarly if K\ = A1A1, K2 = A2A1, we can take ez = C3A1A2, ez = C4A2A2. 
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Figure 7: The splitting of variety into branches under different kinematic configurations. All branches are dimension 
two varieties, except the most degenerate case 

C22 (™,m, m) where 

one branch is dimension four and the other dimension 

three. 



J (0,0) 



7.1 The topology (C32): non-planar crossed box-triangle 

There is only one topology left for type (C), i.e., the crossed box-triangle topology (C32) where three 
propagators containing pure £\, two propagators containing pure £2 and two propagators containing both 
t\ , £2 ■ We use K\ , K% to construct momentum basis e, , and after constructing four linear equations of 
RSPs from seven on-shell cut equations, we can solve, for instance, xi,X2,yi,U2 as functions of four 
ISPs (£3, £4, 2/3, 7/4). The remaining three propagators are then quadratic functions of ISPs. For general 
kinematic configuration, the expression and solution of cut equations are tedious, and we will not explicitly 
write them down here. 

Integrand basis: In general, the variety defined by these three remaining quadratic cut equations is 
irreducible and dimension one. The integrand basis is given by monomials of four ISPs dressed with cor- 
responding propagators. Using Grobner basis under ISPs ordering (j/4, j/3, X4, X3) and the renormalization 
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conditions 

d( Xi )<b, Yl d (^)<4, Y E d (^)<6, (7.1) 

all ISPs of x all ISPs of y all ISPs of x all ISPs of y 

we can get integrand basis in various kinematic configurations. There are all together four kinds of inte- 
grand basis depending on the massless limits of K\,K 3 since we have chosen K\,K 3 as the generators of 
momentum basis. For all kinematic configurations with K\,K 3 massive, the integrand basis contains 38 
elements given by 

&C32 = {!> x 3 , x 2 3 , xl, x 3 , xf , x% x 4 , x 3 x 4 , x§x 4 , X3X4, X3X4, x|x 4 , 2/3, X 3 y 3 , ^32/3, £32/3, 
x 3 yz, x\y 3 , x 4 y 3 , x 3 x 4 y 3 , x 2 3 x 4 y 3 , x\x 4 y 3 , x\x 4 y 3 , y\, x 3 yj, x 4 y\, 2/| , x 3 y|, x 4 yf , 

vhx 3 y 3 ,x 4 y 3 ,y 4 ,x 3 y 4 ,y 3 y 4 ,yly 4 ,yly 4 ) • (7.2) 

For all kinematic configurations with K\ massless while K 3 massive, the integrand basis still contains 38 
elements, and is given by replacing 9 elements in B^ 32 

#C32 = - {^3^4, ^§^4, ^|^4, x\x 4 , x|x 4 , X 3 X 4 2/ 3 , x|x 4 2/ 3 x|x 4 2/ 3 , xjx 4 y 3 } 

+{x 4 , xl, xl, xl, xl, x\y 3 , x 4 2/ 3 , x 4 2/ 3 , x 4 2/ 3 } . (7.3) 

For all kinematic configurations with K\ massive while K 3 massless, the integrand basis contains 38 
elements, and is given by replacing 6 elements in B^ 32 

#C32 = #C32 - {^42/1,^42/1, x 4 yj , y 3 y 4 , y\y±, yly^S + {y\,x 3 y\,y\,x 3 y\,y\, x 3 y\] . (7.4) 

Finally for all kinematic configurations with both K\,K 3 massless, the 38 elements of integrand basis are 
given by replacing fifteen elements in Bq 32 

&C32 = &C32 ~ {^3^4,a:i2;4,X^X 4 ,x|x 4 ,X^X 4 ,X3X 4 2/3,x|x 4 2/3,X^X 4 2/3,x|x 4 2/3, 

x 4 y\, x 4 yl,x 4 y\,y 3 y 4 , y 3 y 4 , y 3 y 4 } + {x\, x\, x\, x 4 , x 4 , x\y 3 , x\y 3 , x\y 3 , 
xly 3 ,y\,x 3 y\,y\,x 3 y\,y\,x 3 y\\ . (7.5) 

To discuss the structure of variety, we again use the notation C32^'p^ where now U, P could either 
be K$,Kq or representing corresponding K^,Kq absent. L will be m if at least one of K\,K.2 being 
massless, and R will be m if K 3 is massless, while N would be m if K 4 is massless, otherwise they will be 
denoted as M. 

The number of branches under various kinematic configurations is summarized in table |]. For each 
kinematic configuration, one should use all branches to fit all 38 coefficients of integrand basis. But we 
can use branch-by-branch polynomial fitting method to simply calculations. 

Variety with one branch: For the most general kinematics C32^^^\ the variety is irreducible 
with dimension one. All 38 coefficients should be fitted at the same time by using this branch. 
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(K 5 /K 6 ,0) 


(0,0) 


(M , M, M) 


1 
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4 


(m, M, M), (M, m, M), (M, M, m) 


2 


4 


6 


(m, m, M), (m, M, m), (M, m, m) 


4 


6 


6 for (M,m,m)/8 


(m, m, m) 


8 


8 


8 



Table 1: Number of branches in various kinematic regions for non-planar crossed box-triangle topology. The 
kinematic configurations are denoted by 032^'^'^. 



Variety with two branches: For kinematic configurations 

p „ 9 (M,M,m) p „ 9 (M,m,M) po ( m .^ M ) r , (M,M,M) ,„ R x 

the variety is given by two branches with dimension one. These branches will intersect at points. More 
explicitly, for C32|^^*^, two branches intersect at two isolated points, while for C32^^jf^\ C32^"^^ 

and 032^^*^, two branches intersect at four points. 

Variety with four branches: For kinematic configurations 

rii „(M,m,m) p„ 9 (m,m,M) p„ 9 (m,M,m) p „ 9 (M,M,m) 

^ 6Z (K 5 ,K 6 ) ' ^ 6Z (K 5 ,K 6 ) ' ^ 6Z (K 5 ,K 6 ) ' ^ 6Z (K 5 /K 6 ,0) ' 

po (M,m,iW) p„ (m,.M-,M) p „ 9 (M,M,M) ,„ „x 

^ JZ (Jf 5 /^6,0) ' ^ JZ (^ B /-fsr6,0) ' ^ Z (0,0) ' { '-'> 

the variety is given by four branches with dimension one. The intersection pattern among four branches 
Vi,V 2 ,V 3 ,V 4 can be shown as follows. For C32^'^^\ C32^^^ and C32^^^s, the only non-zero 

intersections are V 1 nV 2 = V 3 C\V 4 = (0|2), V 1 C\V 3 = V 2 nV 4 = (0|1). For C32[JJ^™ } and C32j^' M) , the 
non-zero intersections are VinV 2 = V 2 C\V 4 = V 4 C\V 3 = V^n Vi = (0|2). For C32^^*L and C32[JJ"^^, 
non-zero intersections are given by V 1 DV 4 = (0|2), Vi DV 2 = V 2 DV 4 = (0|1), V 1 DV 3 = V 3 DV 4 = (0|1). 
Variety with six branches: For kinematic configurations 

p, 9 (m,M,ra) p„ 9 (m,m,M) p „ 9 (M,m,m) p„ 9 (M,M,m) 

^ 6Z (K 5 /K 6 ,0) ' ^^(KrJKe,®) ' ( -" 3Z (E- 5 /^6,0) ' ^ Z (0,0) 

n „ n (M,m,M) POT (m,M,M) nQO (M,m,m) /- q n 

UciZ (0,0) ' UCSZ (0,0) ' UciZ (0,0) ' ^ 

the variety is given by six branches with dimension one. These branches again intersect at points. For 

C32 (SS0) aIld CS2 W&> 6aCh ^ ° f ^ V ^ ( V *> V ^> ^ V3 ^ ^ V5 ^ 

(V 4 ,Vq) intersects at one single point. For 032^^^^, C32^'^ ) ' rn \ C32^'^' M ^ and C32^'™^' m \ each 

pair of (Vi, Vi) and (V 4 , Vi) for i = 2, 3, 5, 6 intersects at one single point. Finally for C32^'^' M \ non-zero 
intersections are given by V 2 D Vi = V 2 n V 4 = V 3 n V 5 = V 3 n V G = (0|2) and Vi DV 5 = V 4 DV 6 = (0|1). 
Variety with eight branches: For kinematic configurations 

p,n(ro,m,m) p„ 9 (m,m,m) p„„(m,M,m) p„„(m,m,M) p„ 9 (m,m,m) ,„ Q x 

^Vs,^) ' ^^(^5/^6,0) ' ^(0,0) ' ^(0,0) ' ^^(0,0) ' ^'- y > 
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(a) 



C32 



(b-1) 



C32 



(b-2) 



C32<™'™;^ } 



( -ig i) (m.,M,M) 
(C-l) 




Figure 8: Intersections of branches for various kinematic configurations of non-planar crossed box-triangle topology 
(C32). Each branch Vi is represented by closed loop and denoted by i, while black dot is the intersecting point. 
Kinematic configurations for each pattern are listed below each diagram. 



the variety is given by eight branches with dimension one. There will be single intersecting point for each 
pair of following ten combinations: (Vi,V 2 ), (Vi.Vg), (V U V 5 ), (V 2 ,V 4 ), (V 3 ,V 4 ), (V 4 ,V 8 ), (V 5 ,V 6 ), (V 5 ,V 7 ), 
(V 6 ,V 8 ) and (V 7 ,V S ). 

The intersection pattern of branches for each kinematic configurations are shown in Figure (^|) . 
7.2 The topology (B41): planar penta-bubble 

There are six propagators for this topology and up to five external legs. From on-shell cut equations of 
these six propagators we can get three linear equations for pure l\, and reduce four RSPs (xi, x 2 , #3, £4) 
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to one. Exception happens when K 4 = K§ = 0, D 4 = (l\ — K\ — K2 — K 3 ) 2 = i\ from momentum 
conservation, and the independent linear equations containing pure i\ reduce to two. In this case we get 
two ISPs from x\. There is no linear equation for pure £2, so all four yi are ISPs. Adding them together 
there will be 5 ISPs (or 6 ISPs for the case K 4 = K 5 = 0). 

We use Ki,K 3 to construct momentum basis ej, and expand all kinematics with this basis. After 
solving linear equations we can express the remaining three quadratic equations with ISPs. We can 
generate Grobner basis from these three quadratic equations, and get integrand basis from all possible 
monomials that satisfying renormalization conditions 

d( Xi )<5, d (^)<2, %*)< 5 - ( 7 - 10 ) 

all ISPs of x all ISPs of y all ISPs of x all ISPs of y 

For instance, if we solve the linear equations of kinematic configurations B41t K4jK& \, B41(^ 4 /^ 5 ) and get 
ISPs (X4, 2/2)2/3) 2/4)) then we can get 18 elements of integrand basis by using Grobner basis under ordering 
(2:4, 2/1)2/2)2/3) 2/4)- These elements will change according to kinematics of K±,K 3 since we have chosen them 
as generators of momentum basis ej. We have three kinds of integrand basis for kinematic configurations 
B41(^ 4 ^ 5 ), B41( K4 / K& y The first kind is suitable for all kinematic configurations of B41(# 4 ^ 5 ), or K 3 
massive while others arbitrary for B41(^ 40 ), or K\ massive while others arbitrary for B41(^ 5j0 ), and these 
18 elements are given by 

#541 = I 1 ) Z4, 2/1, 2/2, 3:42/2, 2/1) 2/3, 2/12/3, 2/22/3, 2/3) 2/4,2?42/4, 2/12/4, 2/22/4, X4V2V4, 2/32/4, 2/i ^42/1} • ( 7 - n ) 

The second kind is suitable for kinematic configurations K\ massless while others arbitrary for B41(^ 6 ), 
and the 18 elements are given by replacing one element in B J B4l 

Bmi = ~ {XAV2VA} + {xavI} ■ (7-12) 

The third kind is suitable for kinematic configurations K 3 massless while others arbitrary for B41(#- 4 ,0), 
and the 18 elements are given by replacing three elements in Bq AI 

b bIi = Bmi ~ {2:42/2, 2/2, 2:42/22/4} + {2:42/1, 2/i, 2:42/1} • (7.13) 

For all kinematic configurations of B41( ) where ISPs are given by six variables (x 3 , x 4 , 2/1,2/2,2/3,2/4), we 
can get 83 elements for integrand basis 

#B4i = I 1 ) 2:3 , x\ , x\ , x\ , x\, x 4 , x\, xl, x\, x 4 , 2/1 , x 3 2/i , x\y x , xjjyi , x\y x , x 4 y x ,x\y u x\y x , x\y x , 
2/i, x 3 yf, xjyl,xlyf, x 4 y\,x\y\, x 4 y 2 ,y 2 , x 3 y 2 , x\y 2 , x\y 2 , X32/2, 2:42/2, 2:42/2, 2:42/2, xjy 2 , 
2/2, 2:32/2, x\y\,x\y\, x 4 y\,x\y\, x 4 y%,y 3 , x 3 y 3 , x\y 3 , x%y 3 , xjy 3 , x 4 y 3 , 2/12/3, 2:32/12/3, 
2:12/12/3, 2:|2/i2/3, 2:42/12/3, 2/22/3, 2:32/22/3, 2:|2/22/3, 2:|2/22/3, x 4 y 2 y 3 , 2/1, x 3 y 3 ,xly 2 , xly 3 ,x 4 y 3 , 
?/4, 2:42/4, 2:42/4, 2:42/4, 2:42/4, 2/12/4, 2:42/12/4, 2:42/12/4, 2:42/12/4, 2/22/4, 2:42/22/4, 2:42/22/4, 2:42/22/4, 
2/32/4, 2:42/32/4,2/4)2:42/4) 2:12/4)2:42/4} • (7-14) 
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The reason we have 83 elements instead of 18 is that, for B41(^ 4 and B41(^ 4 /^ 5 ), 24 is determined by 
quadratic equation, i.e., the maximal power of 24 is two, while for B41( 0j0 ) we can have d(xs) + ^(24) < 5. 

In order to simplify the fitting of coefficients, we need to discuss the branch structure of variety 
defined by the three quadratic equations. For B41( Ki Kr ^ and B41( K4 / Kr> ^, there is a quadratic equation 
of single variable X4, and we can always get two solutions of X4 in C-plane no matter what the momentum 
configuration of K±, K2, K3 is. Thus there will always be two separate branches characterized by two 
solutions 2 4 x ,2 4 2 . For ~B41r K4K& \, these two branches will not split further. Since the three quadratic 
equations can fix 3 ISPs, the two branches can parameterized by two free parameters, and they are 2- 
dimensional. By using each branch we can detect 9 coefficients of integrand basis, and since there is no 
intersection between two branches, we can detect all 9 + 9 = 18 coefficients by using both branches. For 

B41(^ 4 /^ 5)0 ), each branch will split further into two branches, so there will be in total four branches: 

r r 

Vi,V2 characterized by x 4 * and Vs,V4 characterized by 2 4 2 . Each branch can detect 6 coefficients. The 

two branches characterized by 2 4 l will intersect at one-dimensional variety, and we have (1|3) = V\ D V2 
as well as (1|3) = V3 PI V4. So by using two branches of each x 4 * we can detect 6 + 6 — 3 = 9 coefficients, 
and in total 9 + 9 = 18 coefficients by using all 4 branches. 

For kinematic configurations B41/ 0j0 ), the quadratic equation that containing only t\ has two variables 
(23, 24), and we could not divide variety into two separate regions. If all K±, K2,K% are massive, the variety 
is given by two branches. The three quadratic equations can fix three ISPs, and the remaining three can 
be chosen as free parameters. So each branch is 3-dimensional, and the intersection of these two branches 
is 2-dimensional. Each branch can detect 58 coefficients, while the intersection is (2| 33) = V\ D V2. If at 
least one of K\,K% is massive, the variety will split into four 3-dimensional branches Vl,V2 and V?,^. 
By using V\ or V3 we can detect 28 coefficients of integrand basis, while using V2 or V4 we can detect 
36 coefficients. The intersection of these four branches is 1-dimensional, and it can detect 3 coefficients. 
Intersection of every three branches is also the same 1-dimensional variety as intersection of four branches. 
For intersection of every two branches, (Vi,V£) and (V^Vi) are inherited from the intersection of four 
branches, which is 1-dimensional variety. The intersection for (Vi, V4) or (V^Vjj) is 2-dimensional, and we 
have (2 1 6) = V\ D V4, (2|6) = V2 D V3. The intersection for (Vi, V2) or (V3, V4) is also 2-dimensional, and 
we can detect 18 coefficients by using each intersection. 

For the special kinematic configuration of B41( 00 ) with both Ki,K^ massless, the three quadratic 
on-shell cut equations reduce to 

X3X4 = , 2/12/2 + 2/32/4 = 0, 232/4 + 242/3 = 0. (7.15) 

Besides the ordinary four branches 
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there is also another embedded branch given by the ideal 

V 5 : {x A y 3 + x 3 y 4 , 2/12/2 + 2/32/4, xl,x 3 x 4 , £3} 



(7.17) 



Each of the ordinary branches can detect 28 coefficients, while V5 can detect 37 coefficients. These five 
branches intersect at one single point, and intersection of every four branches is also the same point. For 
intersection of every three branches, we have (1|6) = (Vi,V 2 ,Vx,), (1|6) = (V3, 14, U5), ( 1 1 3) = (V"i, 1/3,1/5), 
(1|3) = (V2, 1/4,^/5), and they are all located at different 1-dimensional variety. The other intersection of 
every three branches are inherited from the same point as intersection of five branches. For intersection 
of every two branches, ( Vi , V4) , (V 2 ,V 3 ) are located at the same point as intersection of five branches, 
and (Vi,V 3 ), (1/2,14) are the same 1-dimensional variety as intersection of (V^., V3, V5) and {y 2 ,V A ,V^) 
respectively. Intersections of other combinations are 2-dimensional, and we have (2 1 12) = V\ fl V5, (2| 12) = 
V 2 n V 5 , (2|12) = V 3 n V 5 and (2|12) = V 4 D V 5 , (2|15) = Vi n V 2 , (2|15) = V3 D V4. They are all embedded 
in different 2-dimensional regions. 

7.3 The topology (B33): planar double-box 

This topology has been discussed in details in many other papers [21-24], here we will briefly summarize 
some results. We use Ki,K A to construct momentum basis ei,i = 1,2,3,4 and all kinematics can be 
expanded by this basis. There are seven on-shell cut equations, and among them we can get four linear 
equations and three quadratic equations with multi- variables. We can use four linear equations to reduce 
8 RSPs to 4 ISPs. Since there are two linear equations containing only £1, and two containing only £2, 
we can get, for instance, 4 ISPs (x 2 , X4, yi, 2/4). Then Dq = becomes a conic section of (x 2 ,x A ), Do = 
becomes a conic section of (2/1,2/4) an( i A) = is a multi- variable quadratic equation of (x 2 , £4, 2/1, 2/4). 
These conic sections will degenerate when K\ or K 2 goes into massless limit for Dq = 0, and K 3 or K A 
goes into massless limit for D$ = 0, and variety will split into branches accordingly. 
Under the renormalization conditions 

d( Xl )<4, d (^)< 4 ' E E %0<6, (7.18) 

all ISPs of x all ISPs of y all ISPs of x all ISPs of y 

we can get 32 elements for integrand basis after dividing all possible monomials by Grobner basis that 
generated from three quadratic cut equations with ordering (x 2 , X4, yi, 7/4). The elements under different 
kinematic configurations will be different, and they can be arranged into four kinds according to massless 
limits of K\,K4, which we have chosen as the generators of momentum basis ej. If Ki,K 4 are massive, 
integrand basis is given by following 32 elements 

Bm3 = i 1 . x 2 , £4 , 2? 2 x 4 , x\ , x 2 x\ , x\ , x 2 x\ , x\ , y 1 , x A yi , x\yi ,x\yi, x\y\ ,y±,x 2 y4, x 4 y 4 ,x\y±, 

x\y±, x\y A , yiy A , x A y x y A , yj, x A y\,y x y\,x A y\y\, y\,x A y\,y x y%x A y x y\, yj, x A y\} . (7.19) 
If K\ is massless while K A is massive, integrand basis is obtained by replacing 6 elements in -6g 33 as follows 
#B33 = ^533 - {^2x4, x 2 x\, x 2 x\, x\y x , x\y x , x\y{\ + {x 2 2 , x%, x 2 , x 2 2 y A , x\y±, x 2 y A } . (7.20) 
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If K 4 is massless while K\ is massive, integrand basis is obtained by replacing 6 elements in B J m3 as follows 
b bs3 = B m3 ~ iviVi, x 4 ym, yiyj,x4yiyl,yiyl, x 4 yiyl} + {yf, x 4 yj, y\, x 4 y\, yf, x 4 yf} . (7.21) 

Finally if both K\,K 4 are massless, integrand basis is given by replacing 12 elements in 6^33, which is 
exactly the 6 elements from the second kind of integrand basis plus the other 6 elements from the third 
kind of integrand basis 

nIV r-r 234 234 2 23 34 4 

■^833 = {*■■> x 2i x 2i x 2t x 4 , x 4 , x 4 , x 4 , y±, x 4 y±, Hi , x 4 yi , Hi , x 4 yi , Di , x 4 y± , y 4 , 

X2V4, x\y 4 , x\y 4 , x 2 y 4 , x 4 y 4 , xjy 4 , x 4 y 4 , xjy 4 , yj,x 4 yj,y 4 , x 4 y 4 ,yj,x 4 yj} . (7.22) 

In order to describe branch structure of variety, we use the notation B33^pj, where again U,P 
could be either Kx,,Kq or 0, and L is denoted by m if at least one of Ki, K2 (or K 3 ,K 4 for R) is massless, 
otherwise it will be denoted by M. For general kinematic configuration B33^' A ^.\, the variety is irreducible 
with dimension one. For kinematic configurations B33fc^K, B33^"^, B33^^ g ^ and B33/^'^f\ the 
variety splits into two branches. Each branch can detect 17 coefficients of integrand basis, and two branches 
intersect at two points, which exactly gives 17 + 17 — 2 = 32 coefficients when using both two branches. 

For B33^"^ g N, the variety is given by four branches. Each branch can detect 9 coefficients. There is 
no intersection for four branches or every three branches, while each of following pairs (Vi, V2), (V2,V 3 ), 
(V3, V 4 ) and (V 4 , V±) intersects at a single point. Thus when combining them together we get 9 x 4 — 4 = 32 
coefficients. For kinematic configurations 633^^^ , B33& 

&,e>v B33 (0'0) } and B33 (S'0) } ' the variet y 

is also given by four branches. Branches V±, V3 can individually detect 5 coefficients while branches V2,V 4 
can individually detect 13 coefficients. The non-zero intersections among branches are still single points 
between following pairs (Vi,V 2 ), (V 2 ,V 3 ), (V 3 ,V 4 ), {V 4 ,V{). 

For kinematic configurations B33^^ g ^ and B33|™'™^, the variety is given by six branches. Among 
these six branches, V±,V 4 can detect 9 coefficients while V2,V3,V5,Vf3 can detect 5 coefficients. Non-zero 
intersections exist only for following pairs of two branches (V\, V2), (V2, V3), (V3, V 4 ), (V 4 , V5), (V5, Vq) and 
(Vq, Vi), and each intersection is a single point. So by using all six branches we can detect 2x9+4x5—6 = 32 
coefficients. 

7.4 The topology (B32): planar box-triangle 

This box-triangle topology has six propagators: three propagators containing only £\, two containing only 
£2 and one containing both £\, £ 2 - We can get two linear equations for (a?i, X2, x 3 , x 4 ) and one linear equation 
for (2/1,2/2,2/3,2/4), and reduce 8 RSPs to 5 ISPs. 

We use K\,K 3 to construct momentum basis ej, and after solving the linear equations, for instance, 
we can get a set of ISPs (2:3, x 4 , y±, y 3 , y 4 ). By using all possible monomials of ISPs that satisfying renor- 
malization conditions 

d(xi)<4, %*)<3, E E d ^)^ 5 ( 7 - 23 ) 

all ISPs of x all ISPs of y all ISPs of x all ISPs of y 
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and the Grobner basis that generated from the remaining three quadratic equations under ordering 
(x 3 , x A , yi, y 3 , 2/4), we can get integrand basis for any momentum configurations. For B32(^- 4j ^ 5 ) and 
B32(^- 4 /^ 5j0 ), we can always get 69 elements for integrand basis, but these elements may be different. 
The difference can be classified by momentum configuration of Ki,K 3 , and there are in total 4 kinds of 
integrand basis. The first kind is for all configurations with K A ,K 3 massive, and the 69 elements are given 
by 

Bm2 = { 1 > x 3 , xa , X3X4 , x\ , x 3 x\ , x\ , x 3 x\ , x\ , y 1 , x A yi , x\yi ,x\y±, x\y x ,y 3 ,x 3 y 3 , x A y 3 , x\y 3 , 
x\y 3 ,x\y 3 ,yiy 3 ,x A yiy 3 ,yl, ar 3 y|, x A y\,x\y\,x\y\,yxy\, x 4 yiy%, yf, x 3 y%, x 4 yf, a^yf, 
y 4 , x 3 y A , x A y A , x 3 x A y A , x\y A , x 3 x\y A , x 3 A y A , x 3 x\y A , xjy A , yiy 4 , x A yiy A , x\y x y A , x\y x y A , 
V3V4, x A y 3 y A , x 2 A y 3 y A , x 3 A y 3 y A , y A y 3 y A , x A y A y 3 y A , y 3 y A , x A y 3 y A , x\y\y A , yj,x 3 yj, x A y 2 A , 
xjyl, x\yl,yiyl,x A yiyl, y 3 y\,x A y 3 y\,x\y 3 y\,y\,x 3 y\, x A y\,x\y\\ . (7.24) 

The second kind is for configurations with K\ massless while K 3 massive, and the 69 elements are given 
by replacing 15 elements in B T m2 

&B32 = &B32 ~ {x3X4,x 3 xl,xsxl,xly 1 ,x A y 1 ,xiy 1 ,xlyl,xly 3 ,xlyl,xsX4y4,x 3 xly4, 

3 2 3 22 T,r2 3 4 2 3 4 2232 

x 3 x A y A , x A yiy A , x A y A y A , x A y 3 y A \ + {x 3 , x 3 , x 3 , x 3 y 3 , x 3 y 3 , x 3 y 3 , x 3 y 3 , x 3 y 3 , 

232 3 4 22322 3i orr\ 

x 3 y 3 ,x 3 y A ,x 3 y A ,x 3 y A ,x 3 y A ,x 3 y A ,x 3 y A ) . (7.25) 

The third kind is for configurations with K 3 massless while K\ massive, and the 69 elements are given by 
replacing 12 elements in £>£ 32 

#B32 = ®B32 ~ {V3V4, x A y 3 y A ,xly 3 y A ,xly 3 y A ,y 1 y 3 y A ,x A y 1 y 3 y A ,y 3 2 y A ,x A y 3 y A ,xlyjy A , 
V3V A ,x A y 3 y 2 A ,x 2 A y 3 y 2 A } + {yf,x A yf, yf, x A yf, xjy^, x\y x y 3 , yjy 3 , x A yjy 3 , 
x 2 A yiyl,ylyA,x A yly A ,x 2 A yiy 2 A } . (7.26) 

The last kind of integrand basis is for configurations with K A ,K 3 massless, and the 69 elements are given 
by replacing 23 elements in £>£ 32 

#532 = #532 - {^3^4, x 3 x% x 3 x 3 A , x 2 A y A , x 3 A yx, xjyx, x 3 x A y A , x 3 x 2 A y A , x 3 x 3 A y A , x\y x y A , x\y A y A , 
J/33/4, x A y3VA, xjy 3 y A , x\y 3 y A , y±y 3 y A , x A yiy 3 y A , y 2 y A , x A y\y A , x 2 A y 2 y A , y 3 yj,x A y 3 yj, 

2 2t , r 2 3 4 2 23 32 3 4 2 2 2232 

x A y3V A } + {x 3 ,x 3 , x 3 ,y l ,x A y 1 , y A ,x A y 1 ,x 3 y 3 , x 3 y 3 , x 3 y 3 , y A y 3 , x^y^x^x^, 

232342 222322 3i in q>t\ 

x 3 y 3 ,x 3 y A , x 3 y A ,x 3 y A , y x y A , x A y x y A , x 3 y A ,x 3 y A ,x 3 y A \ . (7.27) 

The integrand basis for B32( 0i0 ) can be distinguished by kinematics of K 3 . If K 3 is massive, we still get 69 
elements, while if K 3 is massless, we can get 77 elements instead of 69. The number of elements changes 
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because in this special momentum configuration, the triangle sub-loop is vanishing triangle, and variety is 
3-dimensional. While in other momentum configurations the variety is 2-dimensional. The 69 elements for 
K 3 massive case are given by 

&B32 = {^^3,xl,xl,xj,x i ,xl,xl,xj,y 1 ,x 3 y 1 ,xly 1 ,xly 1 ,X3y 1 , x A yi ,x\yi, x\y\ , x\y x , y 3 , 

^32/3, x\y 3 , x\y 3 , x 3 y 3 , x 4 y 3 , 2/12/3, x 3 ym, xjym, x 3 2/i2/3, x 4 yiy 3 , y\, x 3 y% x\y\,x\y\, 

xwl, 2/12/3 > x myh ^32/12/3, xmvh vh x zyh x lyh x ±yl,yA, £42/4, x\y A , x\y±, xjy 4 , 

2/12/4, x 4 yW4, x\y x y^ £42/12/4, y 3 y 4 , x 4 y 3 y 4 , 2/12/32/4, a?42/i2/32/4, yfy 4 , x 4 y 3 y 4 , yj,x 4 yj, 

22322 2222 2332 3i oo\ 

x 4 y 4 ,x 4 y 4 ,y 1 y 4 ,x 4 yiy 4 ,x 4 y 1 y 4 ,y 3 y 4: ,X4y 3 y 4 ,y 4 ,X4y4 : ,x 4 y 4 } , (7.28) 

and this kind of integrand basis is different from the previous four kinds of B32^ Ki K ^ and B32(^ 4 /^ s ). 
For configurations with K 3 massless, the 77 elements are given by 

#532 = {hx3,xl,xl,xj,X4,xl,xl,xi,yi,x 3 yi,xlyi,xlyi,xiyi,X4yi,xlyi,xlyi, 

4 2 22232 222323 323 323 

x 4 yi , 2/i , ^32/i ,x 3 y 1: x 3 y x , x 4 y x , x 4 y 1 , x i y 1 , y 1 , x 3 y x , x 3 y l , x 4 y 1 , x 4 y l , y 3 , 
^32/3, x ly 3 , xly 3 ,x^y 3 ,X4y 3 ,y 1 y 3 ,x 3 y 1 y 3 ,xly 1 y 3 ,xly 1 y 3 ,X4y 1 y 3 ,yfy 3 , x 3 yfy 3 , 

22 2 2 22232 2 2223 323 

£32/12/3 , X4yi 2/3, 2/3 , ^32/3 , £32/3 > ^32/3 > 2/12/3 , £32/12/3 > £32/12/3 , 2/3 , £32/3 , £32/3 > 2/4, 

£42/4, Z42/4, x|y 4 , x 4 y 4 , yij/4, x 4 2/i2/4, x\y\y±, x\y\y 4 , yfy^, x 4 yly 4 , x\y\y 4 , y\, 

222322 222332 3i on\ 

^42/4! ^42/4, ^42/4; 2/12/4^42/12/4, ^42/12/4, 2/4, ^42/4, ^42/4} ■ ( 7 - 29 ) 

After getting integrand basis, we need to discuss the branch structure of variety. For the most general 
momentum configuration B32(^- 4 with all external momenta massive, the variety has only one irreducible 
branch, but for some kinematic configurations it will split into many branches. We will use the notation 

IT n\ 

B32^'pj, where as usual U, P could be K4,K§ or 0, while L is denoted by m if at least one of ifi , .K2 (or 
K 3 for R) is massless, otherwise it is denoted by M. The number of branches from the splitting of variety 
can be summarized in table (^). 



(L,R) 

(U,P) 


(K4,K 5 ) 


(K 4 /K 5 ,0) 


(0,0) 


(M, M) 


1 


2 


2 


(m,M), (M,m) 


2 


4 


4 for (m,M); 2+1 for (M,m) 


(m, m) 


4 


6 


2+1 



Table 2: Number of branches in some kinematic regions for planar box-triangle topology (B32). Generally each 
branch is 2-dimensional variety, but for momentum configurations B32^^ and B32|™'™^, there is an extra branch 
of dimension one. We write it as 2 + 1 to emphasize the difference. 

Variety with one branch: For kinematic configuration B32^^K, the variety is irreducible with 
dimension two. All 69 coefficients should be fitted at the same time by using this branch. 
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Variety with two branches: For kinematic configurations 



Rq9 (m,M) R o 9 (M,m) R o 9 (M,M) R oo(M,M) 
a6Z (K 4 ,K 5 ) ' D6Z (K 4 ,K 6 ) ' D6Z (K 4 /K 6 ,0) ' C,iZ (0,0) ' 

the variety has two branches. For B32^^ s, each branch can detect 38 coefficients, and intersection of 
these two branches is 1-dimensional (1|7) = V\ D Vj. For B32^"^, B32|^5j^) p ^ and B32|^ V ^, each 
branches can detect 42 coefficients, and intersection of these two branches is (1|15) = V\ n Vj. 
Variety with four branches: For kinematic configurations 

U o (m,m) R q 9 (m,M) R „ 9 (Af,m) R oo( m > M ) 

the variety has four branches. For B32fe^"^ each branches can individually detect 23 coefficients, while 
for B32^^ e ^ and B32^'^, each V\ or V4 can detect 17 coefficients, and each Vi or V3 can detect 
29 coefficients. The intersections of branches for these kinematic configurations are following. These four 
branches will intersect at one single point, while intersection of every three branches is also the same point. 
For intersection of every two branches, (Vi, V4) and (V2, V3) intersect at the same single point as intersection 
of four branches, and intersections for other pairs are (1|4) = V\ D V3, (1|4) = V2 D V4, (1|8) = V\ D V2, 
(1|8) = V3 fl They are four different 1-dimensional varieties. For B32^^ g0 j, each Vi or V3 can 
detect 10 coefficients, while each V2 or V4 can detect 36 coefficients. There is no intersection for four 
branches, while (Vi,V2, V4) intersect at one single point, and (V^V?, V4) intersect at another single point. 
There is no intersection between (Vi, V3), and intersection of (V2, V4) is 1-dimensional (1|9) = V2 n V4. For 
other intersections of every two branches, we have (1|4) = V\ fl V2, (1|4) = Vi n V4, (1|4) = V2 n V3 and 
(1|4) = V3 n V4. They are different 1-dimensional varieties. 

Variety with six branches: For kinematic configurations 

non(m,m) 

the variety has six branches Vi, V2, V3, V4, V5, Vej. Each Vi or V4 can detect 10 coefficients, each V2 or V5 
can detect 17 coefficients while each V3 or Vq can detect 23 coefficients. There is no intersection among six 
branches and every five branches. The only non-zero intersection of every four branches is (V2, V3, V5, Vq), 
and they intersect at a single point. For intersections of every three branches, (V^V^Vs), (V^V^Ves), 
(V2,V§,Vq) and (Vs,V^,Vq) will intersect at the same point as intersection of (V2, V3, V5, Vq). (Vi,V2,V3) 
will intersect at different single point, and (V4, V5, Ve) will intersect at another different single point. For 
intersections of every two branches, (V2, Vj>) and (V3, Vq) will intersect at the same point as intersection of 

{v 2 ,v 3 ,v 5 ,Vq), while (i|4) = v 1 n v 2 , (i|4) = v 1 n v 3 , (i|4) = v 2 n Vq, (i|4) = v 3 n v 5 , (i|4) = v 4 n v 5 , 

(1|4) = V 4 n Vq and (1|5) = V 2 n V 3 , (1|5) = V 5 n V 6 . They are all different 1-dimensional varieties. 
Variety with 2+1 branches: For kinematic configurations 

B32&$ , B32<£# , (7.30) 
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the integrand basis contains 77 elements, and the three quadratic equations reduce to 

£3^4 = , y 3 y 4 = , x 4 y 3 + x 3 y 4 = . (7.31) 

There will be three branches. Two branches Vi,V>2 are given by x 3 = 0, y 3 = while 1/1,0:4, 2/4 free 
parameters and X4 = 0, 2/4 = while j/i, x 3 , y 3 free parameters. These two branches are 3-dimensional. The 
third branch V 3 is embedded in these two branches, and it is given by the ideal 

V 3 : {x 3 , X3X4, x\, y|, y 3 y 4 , yj, x 3 y 4 + y 3 x 4 } ■ (7.32) 

Geometrically it is just the 1-dimensional variety x 3 = X4 = y 3 = y 4 = with yi free parameter. Although 
the third branch is the intersection of V\ , V2 geometrically, from the point of algebraic geometry, it is an 
independent branch. Each V\ or V2 can detect 39 coefficients, while V 3 can detect 27 coefficients. Since 
geometrically is the intersection of V±,V2, it is clear that intersections of these three branches or every 
two branches are the same 1-dimensional variety, and we have (1|4) = V\ D V 2 , (1 1 14) = V\ D V3 = V2 PI V3, 

(i|4) = y 1 ny 2 nv 3 . 

7.5 The topology (B31): planar box-bubble 

This topology contains a sub-loop of bubble structure. When K 3 = K 4 = 0, there is no difference between 
propagators Dq = £q and D2 = (£q — K\ — K2) 2 because of momentum conservation. This will effectively 
eliminate one on-shell cut equation. For B31(^- 3j ^- 4 ) and ~B31(k 3 /k 4 ,<2)) there are five independent on-shell 
cut equations, and from which we can get two linear equations for (xi, X2, x 3 , x 4 ). By solving these linear 
equations we can reduce 8 RSPs to 6 ISPs. For B31( 0)0 ), we have four independent on-shell cut equations, 
and we can only construct one linear equation for (xi, X2, x 3 , X4). In this case we get 7 ISPs. 

For B31(^ 3i ^- 4 ) and B31^ 3 IKa,<z>) we can use K\,K2 to construct momentum basis e^, while for B31( 0; 0) 
there are only two external legs, we should choose another auxiliary momentum together with one of K±, K2 
to construct momentum basis ej. By expand all momenta with this basis, we can express cut equations 
with 8 RSPs, and by solving the liner equations we can get, for instance, 6 ISPs (x 3 , X4, yi, y2, y 3 , y 4 ) for 
B31(i<r3,K 4 ), B31 (jft:3/jR - 4i0 ), and 7 ISPs (x 2 , x 3 , x 4 , y x , y 2 , y 3 , Va) for B31 (0i0) . The remaining three quadratic 
equations are expressed with these ISPs. By using all monomials of ISPs that satisfying renormalization 
conditions 

d( Xi )<4, %o<2, E d (^)^ 4 ( 7 - 33 ) 

all ISPs of x all ISPs of y all ISPs of x all ISPs of y 

and the Grdbner basis that generated from three quadratic equations under ordering (x 3 , X4, yi, y 2 , y 3 , y 4 ) 
for B31(^ 3; ^ 4 ), B31(^ 3/Ar4)0 ) and (x 2 , x 3 , x 4 , y±, y 2 , y 3 , y±) for B31 (0j0 ), we can get integrand basis for 
all possible momentum configurations. For all possible momentum configurations of B31(^ 3 ^ 4 ) and 
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B31rK 3 /K 4 ,<2)i * ne integrand basis contains 65 elements given by 

#B3i = {!> £3, xl, xl,xj,x 4 , xl, x\, x\, yi,x 3 yi, x\y x , x%yi, x 4 yi, x\y x , x\y x ,y\,x 3 y\, 

x\y\,x A y\,x\y\, y 2 , x 3 y 2 , x\y 2 , x\y 2 ,x A y 2 ,x\y 2 ,x\y 2 ,y\,x 3 y\, x\y\,x A y\, x\y\, 
yz-,x 3 y 3 ,x\y 3 , x^y 3 , x A y 3 , yiy 3 , x 3 yiy 3 , x\yyy 3 , x A yiy 3 , y 2 y 3 , x 3 y 2 y 3 , x\y 2 y 3 , 
x A y 2 y 3 , y 3 , x 3 y 3 ,xly 3 , x A y\,y A , x A y A , x\y A , x\y A , y x y A , x A y x y A , x\y x y A , y 2 y A , 
x A y 2 y A ,xly 2 y A ,y 3 y A ,x A y 3 y A ,yj,x A yl,x 2 A yl} . (7.34) 

For all possible momentum configurations of B31( ), the integrand basis contains 145 elements given by 

#B31 = I 1 ) X 2, X 3 , X 2 X 3 , xl, X 2 x\, x\,X 2 x\,x\, X A , X 2 X A , X 3 X A ,X 2 X 3 X A , X%X A , X 2 x\x A , X3X4, 

xl, x 2 x\, x 3 x\,x 2 x 3 x\, x\x\, x\, x 2 x% x 3 x\, x\, y\,x 3 yx,x\y\, x\y x , x A yi, x 3 x A yi, 
x\x A yx,x\y\,x 3 x\y\,x\y\, y\, x 3 y\,x\y\, x A y\,x 3 x A y\, x\y\,y 2 ,x 2 y 2 ,x 3 y 2 ,x 2 x 3 y 2 , 
x\y 2 , x 2 x\y 2 ,x\y 2 ,x A y 2 ,x 2 x A y 2 , x 3 x A y 2 , x 2 x 3 x A y 2 , x\x A y 2 ,x\y 2 ,x 2 x\y 2 , x 3 x\y 2 , 
x\y 2 , y\,x 3 y\,x\y\, x A y 2 ,x 3 x A y 2 , x\y\, y 3 , x 2 y 3 , x 3 y 3 , x 2 x 3 y 3 , x^y 3 , x 2 x 2 3 y 3 , x\y 3 , 
x A y 3 , x 2 x A y 3 ,x 3 x A y 3 ,x 2 x 3 x A y 3 ,x 2 3 x A y 3 , x\y 3 , x 2 x\y 3 ,x 3 x\y 3 ,x\y 3 , y\y 3 ,x 3 y x y 3 , 
x\y\y 3 , x A ym, x 3 x A y x y 3 ,x\y x y 3 ,y 2 y 3 , x 2 y 2 y 3 ,x 3 y 2 y 3 ,x 2 x 3 y 2 y 3 , x 2 3 y 2 y 3 , x A y 2 y 3 , 
x 2 x A y 2 y 3 , x 3 x A y 2 y 3 ,x\y 2 y 3 ,y\, x 2 y\,x 3 y\,x 2 x 3 y\,x\y\, x A y 3 ,x 2 x A y 3 , x 3 x A y 3 , x\y\, 
y A , x 2 y A , x 3 y A ,x 2 3 y A , x\y A , x A y A , x 2 x A y A , x 3 x A y A , x\x A y A , x\y A , x 2 x\y A ,x 3 x\y A , x\y A , 
ViV4, x 3 yiy A , x 2 3 yiy A , x A y x y A , x 3 x A yiy A , x\y x y A , y 2 y A , x 2 y 2 y A , x 3 y 2 y A , x\y 2 y A , x A y 2 y A , 
x 2 x A y 2 y A , x 3 x A y 2 y A , x\y 2 y A , y 3 y A , x 2 y 3 y A , x 3 y 3 y A , x\y 3 y A , x A y 3 y A , x 2 x A y 3 y A , x 3 x A y 3 y A , 
x 2 A y 3 y A ,y 2 A ,x 2 y 2 A ,x 3 y 2 A ,x A y 2 A ,x 2 x A y 2 A ,x 3 x A y 2 A ,x 2 A y 2 } . (7.35) 

After getting integrand basis, we should analyze branch structure of variety defined by the three 
quadratic equations. 

Branches of B31(^ 3j x 4 ) : The irreducible variety will split into two branches V\, V 2 when at least one 
of K\,K 2 is massless. These two branches are 3-dimensional, and the intersection of them is 2-dimensional. 
Each branch can detect 37 coefficients, while 9 coefficients can be detected by their intersection. So using 
both branches we can detect 37 + 37 — 9 = 65 coefficients. 

Branches of B31(x 3 /K4,0) : There are two 3-dimensional branches if both K\,K 2 are massive. Each 
branch can detect 46 coefficients, and 27 coefficients can be detected by their 2-dimensional intersection. 
When at least one of K\,K 2 is massless, generally the variety will split into four branches Vi,V 2 ,V 3 ,V A . 
Each V\ or V 3 can detect 22 coefficients, while each V 2 or V A can detect 30 coefficients. Intersection of 
all four branches is 1-dimensional, and we have (1|3) = V\ D V 2 PI V 3 D V A . Intersection of every three 
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branches is the same 1-dimensional variety as intersection of four branches. For intersections of every 
two branches, (Vi, V3) and (V2, V4) will intersect at the same 1-dimensional variety as intersection of four 
branches, and all other intersections are 2-dimensional. We have (2|6) = Vi D V4, (2|6) = V 2 D V3 and 
(2| 15) = V\ fl V2, (2| 15) = V3 n V4. They are different 2-dimensional varieties. However, for the special 
momentum configuration B31(#- 3i0 ) where both K 2 ,K 3 are massless, or B31(^ 4i0 ) where both K\,K4 are 
massless, the three quadratic equations reduce to 

X3X4 = , yiV2 + 2/32/4 = , x 4 y 3 + x 3 y 4 = . (7.36) 

There are in total five branches. Four ordinary branches are given by 



Vi : 


V3 


= , yi 


= , 


X3 


= , X4 , 


V2 ■ 


VA ■ 


free parameters , 


V 2 : 


V3 


= , y 2 


= , 


X3 


= , X4 . 


Vl : 


VA ■ 


free parameters , 


V 3 : 


VA 


= , yi 


= , 


x 4 


= , x 3 , 


V2 ■ 


V3 ; 


free parameters , 


V A : 


VA 


= , y 2 


= , 


x 4 


= , x 3 , 


Vi ■ 


V3 ; 


free parameters . 



The fifth branch is given by the ideal 

V 5 : {x A V3 + x 3 Va, VW2 + V3VA, xl, X3X4, X3} . (7.38) 

All these branches are 3-dimensional, and each Vi, V 2 , V3, V4 can detect 22 coefficients while V5 can detect 
37 coefficients. All five branches intersect at one single point. Intersection of every four branches is also the 
same single point. For intersections of every three branches, (1|6) = VinV^nVj), (1 1 6) = VsnV^nVj;, ( 1 1 3) = 
VinV^n V5 and (1|3) = V^flVin V5, and they are different 1-dimensional variety. For intersections of every 
two branches, (Vi, V4) and (V2, V3) are still the same single point, (Vi, V3) is the same 1-dimensional variety 
as intersection of (Vi, V3, V5), and (V2, V4) is the same 1-dimensional variety as intersection of (V2, V4, V5). 
The intersections of all other pairs are 2-dimensional, and we have (2|12) = Vi fl V 2 , (2| 12) = V3 fl V4 and 
(2|12) = Vi n V 5 , (2|12) = V 2 n V5, (2|12) = V 3 n V 5 , (2|12) = V4 n V 5 . They are all different 2-dimensional 
varieties. By using these five branches, we can detect 65 coefficients of integrand basis. 

Branches of B31( 0i0 ): There are only two external legs, and none of them can be massless, so 
we have only one momentum configuration with both K±,K 2 massive. The integrand basis contains 145 
elements. There are two branches of dimension four, and 110 coefficients can be detected by each of them. 
Intersection of these 2 branches is 3-dimensional, and using it we can detect 75 coefficients. So all 145 
coefficients can be detected by using these two branches. 

7.6 The topology (B22): planar double-triangle 

For the double-triangle topology (B22), we can use K±,K 2 to construct momentum basis. When K3 = 
K4 = 0, K\ and K 2 are not independent, and we use K\ and another auxiliary momentum to construct 
momentum basis. There are five propagators and by using two linear equations Dq — D\ = 0, Dq — D\ = 0, 
we can solve x\,y 2 . So there are six ISPs (x 2l X3, X4, y±, j/3, y 4 ) and three quadratic equations left. 
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This topology has Z2 symmetry between K\,Ki and Z<i symmetry between ^3,^4, we will take the 

IT r>\ 

notation B22^' p j where U, P could be K 3 , K 4 or 0, and L is denoted by m if K\ (or K2 for R) is massless, 
otherwise it is denoted by M. It is worth to notice when K 3 = K A = 0, we have K\ = —K2, thus to get 
non-zero contribution, K±,K2 should be massive. In other words, we do not need to consider kinematic 
configurations B22™'^, B22^'™ and B22™'™. 

By using Grobner basis defined from three quadratic equations under ordering (x2, X3, x 4 , yi, y 3 , 2/4) 
and all possible monomials satisfying renormalization conditions 

Y, d (^)<3, E d( yi )<3, Y d ( x *)+ < 4 , (7.39) 

all ISPs of x all ISPs of y all ISPs of x all ISPs of y 

we can get 111 elements for integrand basis. The explicit form of these elements depends on the kinematics 
of Ki,K 2 , which we have chosen as generators of momentum basis e^. There are in total four kinds of 
integrand basis. For momentum configurations with K\,Ki massive, the 111 elements are given by 

#B22 = {!) x 2, x 3 , x 2 x 3 , x§, x 2 x§, X3, x 4 , x 2 x 4 , X3X4, x 2 x 3 x 4 , X3X4, X4, X2X4, X3X4, xl,yi, 
x 3 y\, xjyi, x\y u x 4 yi, x 3 x 4 yi, x|x 4 yi, x\y u x 3 x\yi, x\yi, y 3 , x 2 y 3 , x 3 y 3 , x 2 x 3 y 3 , 
xjy3, x 2 xly 3 ,xly 3 ,x 4 y 3 ,x 2 x 4 y 3 , x 3 x A y 3 , x 2 x 3 x A y 3 , x 2 3 x 4 y 3 ,x 2 4 y 3 ,x 2 x 2 A y 3 , x 3 x\y 3 , 
x\y 3 , yiy 3 , x 3 yiy 3 , x\yiy 3 , x 4 yiy 3 , x 3 x A yiy 3 ,x 2 4 yiy 3 ,yl, x 2 y%, x 3 y 3 ,x 2 x 3 y 2 ,xly 3 , 
X4y 3 , x 3 x 4 y 3 , xjy 3 ,yiy 3 , x 3 y 1 y 2 ,x 4 yiy 3 ,yi, x 2y 3 i ,x 3 y 3 i , x 4 y%,y 4 , x 2 y 4 , x 3 y 4 , xjy 4 , 
x\y 4 , x A y 4 , X2X 4 y 4 , x 3 x 4 y 4 , xjx A y A , x\y 4 , x 2 x\y 4 , x 3 x\y A , x\y A , yiy 4 , x 3 yiy 4 , x\y\y A , 
x 4 yiy 4 , x 3 x A yiy A , x\yiy 4 , y 3 y 4 , x 3 y 3 y 4 , x%y 3 y A , x 4 y 3 y 4 , x 3 x 4 y 3 y 4 , xjy 3 y 4 , yiy 3 y 4 , 
x 3 yiy 3 y 4 , x 4 yiy 3 y 4 , y 2 y A , x 3 y 3 y 4 , x A y 2 y 4 , yj,x 2 yl, x 3 yj,x 4 yj, x 2 x 4 yl,x 3 x 4 yj,xlyj, 
yiy 4 , x 3 yiy 2 4 ,x 4 yiyl, y 3 yj, x 3 y 3 yj,x 4 y 3 yj, y 4 ,x 2 y 4 , x 3 y 4 , x A y\} . (7.40) 

The second kind of integrand basis is for configurations with K\ massless and K2 massive, and the 111 
elements are given by replacing 19 elements in B I m2 

B 1 ^ = Bm2 ~ {x3X 4 ,X2X 3 x 4 ,xlx 4 ,x 3 xl,x 3 x 4 y 1 ,xlx 4 y 1 ,x 3 xly 1 ,x 3 x 4 y 3 ,x 2 x 3 x 4 y 3 , 
xjx 4 y 3 , x 3 x 2 4 y 3 ,x 3 x 4 y!y 3 , x 3 x A y 2 , x 3 x A y 4 , x\x A y A , x 3 x\y 4 , x 3 x 4 yiy 4 , x 3 x 4 y 3 y 4 , 
x 3 x 4 yj} + {x 2 ,, X2, X2X3, X2X4, x\y 3 , x\y 3 , xlx 3 y 3 ,xjx 4 y 3 ,xjy 3 , x 2 x A y 2 , x\y 4 , 
x%y 4 , X2X3J/4, x 2 2 x 3 y A , x 2 x 2 y 4 , x 2 2 x A y A , x\y\, x 2 x 3 y\,x\y 2 A ) . (7.41) 

The third kind of integrand basis is for configurations with K2 massless and K\ massive, and 111 elements 
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are given by replacing 15 elements in B B22 

Bm2 = &B22 - {y3y4,x 3 y 3 y4,xly 3 y4,x4y3y4,x 3 x 4 y3y4,xly3y4,yiy3y4,x 3 yiy3y4,x l iyiy3y4, 

ylv4, x 3 yjy A , x A y 3 y A , y 3 yj,x 3 y 3 yj, x A y 3 y 2 A } + {yj, x 3 y\,x\y\, x A yf,x 3 x A yj, x\y\, 

Vi,x 3 yf, x A yl,yjy 3 ,x 3 yjy 3 , x A y\y 3 ,y\y A , x 3 yjy A , x A yfy A } . (7.42) 

Finally the fourth kind of integrand basis is for configurations with K±,K 2 massless, and 111 elements are 
given by replacing 33 elements in B I B22 

Bm2 = Bb22 ~ {x3X A ,x 2 x 3 x A ,xlx A ,x 3 xl,x 3 x A yi,xlx A yi,x 3 xlyi,x 3 x A y 3 ,x 2 x 3 x A y 3 ,xlx A y 3 , 
x 3 xly 3 ,x 3 x A yiy 3 , x 3 x A y 3 , x 3 x A y A , x\x A y A , x 3 x\y A , x 3 x A yiy A , y 3 y A , x 3 y 3 y A , xjy 3 y A , 
x A y 3 y A , x 3 x A y 3 y A , x\y 3 y A , yiy 3 y A , x 3 yiy 3 y A , x A yiy 3 y A , yjy A , x 3 yjy A , x A y 3 y A , x 3 x A yj, 



2 2 2i , r 2 3 2 2 2 222 2223 3 3 

V3V A , x 3 y 3 y A , x A y 3 y 4 } + {x 2 ,x% , x 2 x 3 , x 2 x A 7 y 1 ,x 3 y l , x 3 y t ,x A y 4 , x 4 y± ,yi,x 3 y'{, x A yf , 

x\y 3 ,x\y 3 ,x\x 3 y 3 , x\x A y 3 , y\y 3 ,x 3 y\y 3 , x A y\y 3 ,x\y\, x 2 x A y 3 ,x^y A , x\y A , x 2 x 3 y A , 

x\x 3 y A , x 2 x\y A , x\x A y A , y\y A , x 3 yjy A , x A yfy A , x\y\, x 2 x 3 y\, xjyj} . (7.43) 



After getting integrand basis, we continue to present the branch structure of variety. The number of 
branches under different kinematic configurations is summarized in table ||[ 



^\(A#) 

(U,P) 


(K 3 ,K A ) 


(K 3 /K A ,0) 


(0,0) 


(M, M) 


1 


2 


2 


(m,M), (M, m) 


2 


4 




(m, m) 


4 


6 





Table 3: Number of branches in some kinematic regions for planar double-triangle topology (B22). Each branch is 
3-dimensional variety. 

Variety with one branch: For general kinematic configuration B22^^^, the variety is irreducible 
with dimension three. All 111 coefficients of integrand basis can be detected by this branch. 
Variety with two branches: For kinematic configurations 



B22 



(m,M) 
(K 3 ,K4,) 



B22 



(M,m) 
(K 3 ,K 4 ) ' 



B22 



(M,M) 
(K 3 /K 4 ,&) ' 



B22 



(M,M) 

(0,0) 



(7.44) 



the variety is given by two branches V\,V 2 . For B22|^H^y and B22 
efficients, and their intersection is 2-dimensional variety which can detect 31 coefficients. For B22 



,(M,m) 



and B22 



(M,M) 



each branch can detect 71 co- 



(M,M) 
(K 4 /K s ,0) 



, s , each branch can detect 77 coefficients, and their intersection is V\ D V2 = (2|43) . 
Variety with four branches: For kinematic configurations 



B22 



(m,m) 
(K 3 ,K 4 ) ' 



B22 



(m,M) 

(K 3 /K 4 ,0) ' 



B22 



(M,m) 

(K 3 /K 4 ,0) ' 



(7.45) 
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(«) (b) 

Figure 9: Venn diagrams for intersections of four branches, where each ellipse represents one branch. Venn diagram 
(a) is for B22^"^, and Venn diagram (b) is for B22&^- 4 s, B22^*^- 4 y Number means that there is no 
intersection and (d\n) represents that the dimension of that region is d and the number of coefficients detected in 
that region is n. 

the variety is given by four branches. Intersections of branches is expressed by Figure [| 
Variety with six branches: For kinematic configuration 

B22^ 4)0) • (^6) 

the variety is given by six branches Vx,V 2 , Vj, V4, V5, Vg. Among these six branches, four Vi,i = 1,2,3,4 
can detect 29 coefficients and the other two V5, Vq can detect 45 coefficients. All six branches intersect at a 
single point, and intersection of every five branches is also the same single point. For intersections of every 
four branches, most of them are the same single point inherit from intersections of every five branches 
except for the following two combinations of branches (Vi, V3, V5, Vq) and (V 2 , V4, V5, Vq), which intersect 
at 1-dimensional variety, and can be used to detect 4 coefficients. For intersections of every three branches, 
besides the ones that are inherited from intersections of four branches, there are also four pairs (Vi, V 2 , V5), 
(Vx, V 2 , Vq) and (V3, V4, V5), (V3, V4, Vq) that intersect at one-dimensional variety and can be used to detect 
4 coefficients. Intersections of every two branches could be 2-dimensional, 1-dimensional, or single point, 
and they are summarized as 



(2|13) 


= v x n v 2 


= Vxnv 5 


= v 2 nv 6 


= v 3 n y 4 


= v 3 nv G 


= v 4 n v 5 , 




(2|10) 


= v 2 nv 5 


= Vxnv e 


= Vi n v 6 


= v 3 n v 5 








(1|4) 


= Vx n v 3 


= v 2 nv 4 


, (0|1) = 


■- Vx n Vi = 


v 2 nv 3 , 


(h\7) = v 5 nv 6 . 


(7.47) 



One interesting point is that intersection of (V^Vf,) is two 1-dimensional varieties, and to emphasize this 
subtlety we have used (I2) notation. 
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7.7 The topology (B21): planar triangle-bubble 

For this topology we could have two cases H21^ 2 ,K 3 ) an d ^^(k 2 /k- s ,0)- For B21^ 2 x s ), since there are 
three external momenta, we can use K\,K 3 to construct momentum basis. For B21t K2 / K3j(Z) \, there are 
only two external legs, and only one of them is independent. So we need another auxiliary momentum 
together with K\ to construct momentum basis. We do not consider B21( ) since it requires K\ = 0, 
which is tadpole-like structure. For B21(^ 2 /^ 30 ) we also assume that both two external momenta are 
massive for non-vanishing result and do not consider the kinematic configurations where any momentum 
is massless. 

There are four propagators, and from them we can construct one linear equation for X{, and reduce 8 
RSPs to 7 ISPs. For instance, we can solve x\, and express all remaining three propagators as quadratic 
polynomials of seven ISPs {x2,x 3 ,X4,y\,y2,y 3 ,y4)- Using Grobner basis that generated from the three 
quadratic equations with ordering (yi,X2,y2,x 3 , y 3 ,x A , y A ) and all possible monomials of ISPs that satisfying 
renormalization conditions 

all ISPs of x all ISPs of y all ISPs of x all ISPs of y 

we can get integrand basis for different kinematic configurations. The elements of integrand basis depends 
on the kinematics of K\. For B21^ 2 Kz ^ and B2\^ K2 j Kz ^, if K\ is massive, the integrand basis contains 
80 elements given by 

#B21 = {!> X 2 ,X 3 , X2X3, xl, X2X3, X3, X 4 , X 2 X 4 , X3X4, X2X3X4, X3X4, xl, X2X4, X3X4, X4, y\, 

xzy\,x\y x ,X4y\,x\y\,y\, x 3 y\ , x 4 yf, 2/2, £22/2, x 3 y 2 , x 2 x 3 2/ 2 , x§2/2, a?42/2, a?2^42/2, 

x 3 x 4 y2, x\y 2 , y\, x 3 yj,x 4 yj, y 3 , x 2 y 3 , x 3 y 3 , x 2 x 3 y 3 ,xly 3 ,X4y 3 ,x 2 x A y 3 , x 3 x A y 3 , 

xjy 3 ,yw 3 , x 3 yw 3 , x 4 yw 3 ,y2y3, x 2 y2V3, x 3 y 2 y 3 , X4y2y3,yh X2y 3 2 ,x 3 y 3 ,x4y 2 ., y 4 , 

X2V4, x 3 y A , x 2 x 3 y 4 , x\y4, x A y A , x 2 x A y A , x 3 x A y A , x\y A , ym, x 3 yiy4, x A y\y A , 2/22/4, 

^22/22/4, x 3 y 2 y A , X42/22/4, 2/32/4, x 2 y 3 y A , x 3 y 3 y 4 , x A y 3 y A , y\, x 2 y\,x 3 y\, x A yj} . (7.49) 

If K\ is massless for B21(^- 2 integrand basis is given by replacing 8 elements from £>£ 2 i 

#B2i = ~ {x 3 x A ,X2X 3 x A ,xlx4,x 3 xl,x 3 x A y2,x 3 x A y 3 ,x 3 x A y4,X2 2/32/4} 

+{x 2 2 , x\, x\x 3 ,x\x4,x\y 2 , x 2 y\,x\y 3 ,x\y A \ . (7.50) 

The variety defined by the three quadratic equations is irreducible with dimension four for B21(^ 2 ^ 3 ) 
with K\ massive. If K\ is massless, the variety will split into two branches, and each branch can detect 
54 coefficients. Intersection of these two branches is an irreducible 3-dimensional variety, and it can detect 
28 coefficients. So by using both two branches, we can detect 54 + 54 — 28 = 80 coefficients as it should 
be. For B21(k 2 /k 3 ,<z>), K\ should be massive, and the variety has two branches. Each branch can detect 64 
coefficients, and intersection of these two branches is an irreducible 3-dimensional variety, which can detect 
48 coefficients. By using these two branches we can detect 64 + 64 — 48 = 80 coefficients of integrand basis. 
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7.8 The topology (Bll): planar sun-set 

For this topology, since K\ = —K 2 , we use K\ and another auxiliary momentum to construct momentum 
basis. The only possible kinematic configuration is both K\, K 2 massive. There are only three propagators 
and we can not construct linear equation from on-shell cut equations, thus there are 8 ISPs. The three 
quadratic equations can be expressed as 

D = x x x 2 + 23x4 , D = 2/12/2 + 2/32/4 , 

Do = x 2 yi + xiy 2 + x A y 3 + x 3 y 4 + {x 1 + yi)a u + {x 2 + 2/2)012 + an«i2 • (7.51) 

Using Grobner basis with ordering (1/4, y 3 , x 4 , x 3 , y 2 , yi, x 2 , x\) and all possible monomials of ISPs under 
renormalization conditions 

Y, d( Xi )<2, Y «^ 2 ' E E d (Vi)^ 2 ' ( 7 - 52 ) 

all ISPs of x all ISPs of y all ISPs of x all ISPs of y 

we can get 42 elements for integrand basis as 

Bbu = {1, x 1 ,xl,x 2 ,x 1 x 2 ,xl,x 3 ,x 1 x 3 , x 2 x 3 , x\, x 4 , xix 4 , x 2 x 4 , xl,yi,x 1 y 1 ,x 2 yi, 
X3yi,X4yi,yi,y2,xiy 2 ,x 2 y 2 ,x 3 y 2 , X42/2, 2/12/2, 2/1, 2/3, xiy 3 , x 2 y 3 , x 3 y 3 , X42/3, 
2/12/3, 2/22/3, y 3 , 2/4, X12/4, X22/4, a; 4 2/4, 2/12/4, 2/22/4, 2/1} • ( 7 - 53 ) 
The variety defined by the three quadratic equations is irreducible with dimension five. 

8. Conclusion 

In this paper, we use the new technique developed in [27-29] to classify the two-loop integrand basis in 
pure four dimension space-time. Although there are only small number of topologies for planar and non- 
planar two-loop diagrams, the diverse external momentum configurations greatly increase the number of 
integrand basis that we need to discuss. Because the integrand basis and branch structure of variety will 
depend on the topology as well as external kinematics, it is necessary to classify possible sets of integrand 
basis and study the evolution of variety under various kinematic limits. 

The algebraic geometry method, such as Grobner basis method and polynomial division, plays a crucial 
role in our discussion. Using these methods, we are able to present explicit form of integrand basis as well 
as detailed study of varieties, such as branch structures and their intersections. The same method also 
allows us to determine coefficients of integrand basis. 

We must emphasize that our result is only a small step towards the practical evaluation of general 
two-loop amplitudes. As we have mentioned in the introduction, the number of two-loop integrand basis is 
much more than the number of two-loop integral basis and it is highly desirable to reduce integrand basis 
further. One way to do so is to use the IBP-method [18]. However, with the time consuming, it is not 
feasible at this moment. 
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Although we focus on two-loop diagrams in pure four-dimension space-time, the same analysis can be 
applied to(4 — 2e)-dimension for complete answer, or three and higher loop amplitude as demonstrated 
in [29]. It is also an interesting problem to apply these general analysis to real processes. 



Acknowledgement 

We would like to thank Simon Badger, Hjalte Frellesvig, Yang Zhang for many explanations and discussions 
and early participant of this project. B.F would like to thank the hospitality of Niels Bohr International 
Academy and Discovery Center. This work is supported, in part, by fund from Qiu-Shi and Chinese NSF 
funding under contract No.11031005, No.11135006, No. 11125523. 



A. Some mathematical backgrounds 

In this section, we present several mathematical facts that may be useful for determining branch structure 
of non-linear on-shell cut equations. First let us consider the quadratic equation of two variables defined 
by equation 

Ax 2 + Bxy + Cy 2 + Dx + Ey + F = , (A.l) 

with A, B, C not all zero. This equation is usually called conic section. In general, the conic is an irreducible 
one-dimensional variety, however, when the determinant A of following 3x3 matrix 

/ A B/2 D/2\ 



-(B 2 - 4AC)F + BDE - CD 2 - AE 2 



(A.2) 



A = dct B/2 C E/2 
\D/2 E/2 F J 

is zero, the conic splits to two branches. 

Next let us consider the roots of a polynomial 

n 

= (A.3) 

One can use discriminant to determine whether if this polynomial has repeated roots or not. If discriminant 
equals to zero, then there is repeated roots. The simplest example is quadratic equation CL2Z 2 + a\z + ao = 
whose discriminant is D = a 2 — 4aoa2. If D = 0, then this equation has double roots, and the polynomial 
can be written as a perfect square of one factor. The discriminant of n = 2, 3 can be found in many other 
references, and using it we can tell the properties of roots just from coefficients of variable. A special 
interesting example is the quartic function 

f(z) = Az 4 + Bz 3 + Cz 2 + Dz + E , (A.4) 

which can be factorized as 

f(z) ~ (z - z i+ ^){z - z {+ '^)(z - - *(-+>) , (A.5) 
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where z ±,zL are four roots of quartic equation. We want to know if it can be expressed as perfect square 
terms such as 

f(z) = {az + b) 2 {cz + d) 2 . (A.6) 

In other words, we want to know if there are repeated roots or not. The quartic equation can be solved 
by defining 

W 2 C n _ B?> BC D 3B 4 CB 2 BD E 

Then when B = 0, the quartic equation has following solution 



Z = ~4A ±s i 2 ' (A - 8) 

where ± s and ±j can take plus and minus sign independently. If the coefficients further satisfy A 2 — AC = 0, 
then x( s ' + ) = x^ s ~\ and f(x) can be expressed as products of two perfect squares. 

Using above results, we can check whether variety defined by following equations is reducible 

Ax 2 + Bxy + Cy 2 + Dx + Ey + F = , (A.9) 
a{r)x + b{ T )y + c(r) = , (A. 10) 

where a{r) ,b{r) , c{t) are linear functions of (possible free parameter) r. After solving the linear equation 
of x, y and substituting the result into quadratic equation we get 

(a 2 C - abB + b 2 A)y 2 + (a 2 E - abD - acB + 2bcA)y + {a 2 F - acD + c 2 A) = , (A.ll) 

where the coefficients A'(t) = (a 2 C - abB + b 2 A), B'(t) = (a 2 E - abD - acB + 2bcA) and C"(r) = 
(a 2 F — acD + cPA) are now quadratic functions of r. The solution of y is given by 



-B\r) ± y/B'(ry - A A'{t)CJ{t) ^ 



2A'(r) 

and y is a rational function of r when and only when terms inside the square root is perfect square, i.e., 
the quartic function (B'{t) 2 — AA'(t)C'(t)) of r is a perfect square. 
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